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ABSTRACT 


This  report  presents  the  findings  of  a  number  of  individual  investi¬ 
gations  relevant  to  the  physics  of  atomic  and  molecular  processes  in 
the  upper  atmosphere  as  pertinent  to  the  problem  of  the  rate  of  reduc¬ 
tion  of  ionization  in  a  disturbed  atmosphere. 

Calculations  have  been  performed  of  the  rates  for  three-body  ion- 
electron  recombination 

X++e+M-X  +  M 

for  a  variety  of  ions  X+  and  third  bodies  M.  Emphasis  has  been  placed 
on  systems  where  the  third  body  involves  a  molecular  species.  It  has 
been  showr  that  if  the  third  body  has  low  energy  modes  of  internal 
excitation,  large  recombination  rates  may  result.  If  M  is  a  polar  mole¬ 
cule,  very  large  recombination  rates  will  result.  The  calculations  show 
that  the  rate  for  three-body  ion-electron  recombination  is  very  sensi¬ 
tive  to  the  concentration  of  water  vapor. 

Hie  process  of  associative  detachment 

X~  +  Y  —  XY  +  e 

is  an  important  mechanism  for  the  removal  of  negative  atomic  and  molec¬ 
ular  ions.  Only  sparse  accurate  experimental  data  exist,  and  quant' es¬ 
tiva  theoretical  data  exist  only  for  H  +  H  collisions.  Quant itat  .ve 
prediction  of  the  rate  coefficient  for  associative  detachment  rec lires 
accurate  potential  energy  curves  for  the  negative  molecular  ion  XY  and 
the  neutral  molecule  XY,  especially  in  any  regions  of  curve  cr  -s sings 
involving  the  lower  levels  of  the  molecules.  Also  required  --  ;e  certain 
non-radiative  coupling  matrix  elements  between  the  electron  c  states  of 
the  molecules.  The  theoretical  calculation  of  the  necessa-  y  accurate 
potential  energy  curves  and  matrix  elements  represents  a  •  ?latively 
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difficult  theoretical  problem  for  species  of  interest  in  the  atmosphere. 
From  knowledge  of  the  potential  energy  curves,  useful  upper  limits  for 
the  rate  of  associative  detachment  can  be  obtained.  Improved  calcu¬ 
lations  of  potential  energy  curves  Oj,  relative  to  .theoretical  0^  curves, 
have  been  performed  for  the  important  atmospheric  process 

O'  +  0  -  02  +  e 

2 

the  cross  section  could  be  as  large  as  4.6  «aQ  . 

The  cross  section  for  rotational  excitation  in  electron  polar  mol¬ 
ecule  collisions, 

e  +  M(J)  -  e  +  M  (J  +  AJ) , 

where  M(J)  is  a  polar  molecule  in  rotational  state  J,  is  large.  Rota¬ 
tional  excitation  processes  represent  an  important  energy  loss  mechanism 
for  the  thermal iz at ion  of  supra-thermal  electrons  produced  by  an  atmo¬ 
spheric  disturbance.  Since  most  recombination  processes  proceed  much 
more  rapidly  for  low  energy  electrons  than  for  supra-thermal  electrons, 
the  rotational  excitation  process  plays  an  important  role  in  the  rate 
of  reduction  of  ionization  in  a  disturbed  atmosphere.  The  cross  section 
for  this  process  constitutes  an  essential  input  data  point  for  the  studies 
of  three-body  ion-electron  recombination  in  polar  gases.  Detailed  cal¬ 
culations  have  been  performed  for  the  cross  section  for  rotational  exci¬ 
tation  of  carbon  monoxide  by  electron  impact  using  close-coupling  methods. 
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FOREWORD 


This  report  was  prepared  by  GCA  Corporation,  Bedford,  Massachusetts 
under  Contract  No.  DASA01-70-C-0091.  The  research  was  performed  under 
NWER  Subtask  HD028.  Professors  D.  R.  Bates  and  J.  C.  Browne  made  con¬ 
siderable  contributions  to  the  research  program. 
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I.  INTRODUCTION 


A.  detailed  understanding  of  the  behavior  of  the  atmosphere  follow¬ 
ing  a  severe  disturbance  requires  detailed  knowledge  of  cross  sections 
and  rate  coefficients  for  a  highly  diverse  ensemble  of  atomic,  mclecular, 
and  chemical  processes  involving  ambient  atmospheric  species.  Under 
DNA  sponsorship  for  the  past  several  years,  GCA  Technology  Division  has 
been  engaged  in  theoretical  studies  oi  various  selected  atomic  and 
molecular  processes  which  are  relevant  to  the  understanding  of  a  dis¬ 
turbed  atmosphere.  Most  of  the  processes  which  have  been  studied  pertain 
to  the  problem  of  the  rate  of  reduction  of  ionization  following  a  severe 
disturbance  since  this  constitutes  an  important  operational  factor  of 
interest  to  the  sponsoring  agency.  Three  major  criteria  were  employed 
in  the  selection  of  processes  for  detailed  investigation,  as  character¬ 
ized  by  the  following  specific  questions: 

(1)  What  is  the  over-all  importance  or  potential  importance  of  the 
process  in  the  understanding  of  the  de-ionization  of  the  disturbed 
atmosphere:  (In  certain  instances  the  pertinent  rate  coefficients  were 
sufficiently  uncertain  as  to  require  preliminary  investigation  simply  to 
assess  their  relative  i- jportance . ) 

(2)  Can  quantitative  theoretical  calculations  be  performed  to  use¬ 
ful  accuracy  for  atmospheric  species? 

(3)  Are  theoretical  studies  required  either  to  assess  the  accuracy 
or  in  other  useful  ways  to  supplement  available  experimental  data? 

Several  processes  were  evaluated  briefly  and  were  found  not  to  aatisfy 
one  or  more  of  the  above  criteria  and  consequently  were  removed  from 
further  consideration.  Of  the  remaining  candidates,  four  specific  pro¬ 
cesses  were  selected  for  further  detailed  investigation  as  reported 
herein.  Although  each  process  is  believed  to  have  relevance  to  the  gen¬ 
eral  ambient  ionization  reduction  problem  in  conformity  with  item  (1) 
above,  they  are  sufficiently  independent  to  warrant  individual  attention. 
Specifically  the  de-ionization  mechanisms  considered  include:  (1)  three- 
body  recombination,  (2)  ion- ion  neutralization  and  (3)  associative  detach- 
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i  ment.  Finally  some  attention  has  been  given  to  the  theoretical  problem 

ur 

|  of  using  effective  potentials  to  calculate  molecular  wave  functions  and 

i  potential  energy  curves  for  diatomic  molecules.  These  studies  may  then 

l  be  related  to  pertinent  atmospheric  change  transfer  and  scattering  pro- 

l  cesses.  As  a  consequence,  the  detailed  results  of  these  investigations 

I  are  presented  in  Sections  II  through  V  of  this  report  as  discussed  below. 

|  In  Section  II,  three  body  ion-electron  recombination  processes 

s 

b 

|  such  as 

c 

[  X+  +  e+  M-  X  +  K 


are  discussed.  This  represents  an  important  processes  in  the  D-region 
at  its  relatively  low  temperature  and  high  neutral  particle  density, 
especially  for  those  ions  which  do  not  undergo  exothermic  chemical  reac¬ 
tions.  Preliminary  calculations  were  performed  for  various  ions  in  atomic 
third  body  species,  where  proper  account  was  taken  of  the  velocity  de¬ 
pendence  of  the  electron  third  body  scattering  cross  section  (Ref.  1-1). 
Since  most  atomic  species  have  no  low  energy  modes  of  internal  excita¬ 
tion,  only  elastic  electron  third  body  scattering  contributes  and  the 
resultant  recombination  coefficient  is  not  very  large,  although  it  can 
exceed  pure  radiative  recombination.  Atomic  oxygen  represents  an  excep¬ 
tion,  since  fine  structure  transitions  in  the  ground  state  occur  at  low 
energies.  However,  in  the  D-region  the  atomic  oxygen  concentration  is 
not  large  and  molecules  should  constitute  the  major  class  of  third  bodies. 


Calculations  have  been  performed  previously  for  several  homonuclear 
third  bodies  (Ref.  1-2).  Because  of  low  energy  rotational  and  vibrational 
modes  of  internal  excitation,  the  recombination  rate  coefficients  for 
molecular  third  bodies  were  found  to  be  much  larger  than  those  obtained 
for  atomic  third  body  systems.  The  present  calculations  were  performed 
for  polar  molecule  third  body  systems.  Owing  to  the  relatively  long 
range  nature  of  the  electron-permanent  dipole  interaction,  the  cross  sec¬ 
tion  for  rotational  excitation  of  polar  molecules  is  very  large,  so  that 
correspondingly  high  three  body  ion  electron  recombination  coefficients 
result.  The  inelastic  cross  section  for  these  systems  exceeds  the  elastic 
cross  section. 
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These  calculations  demonstrate  that,  since  three  body  ion-electron 
recombination  coefficients  involving  polar  molecule  third  bodies  are 
„.'cen  several  orders  of  magnitude  faster  than  coefficients  for  atomic 
species,  the  over-all  recombination  rate  in  dense  neutral  gases  is  a 
sensitive  function  of  the  polar  molecule  concentration.  The  theoretical 
studies  conducted  to  date  has  provided  sufficient  information  to  esti¬ 
mate  the  recombination  coefficient  for  an  extensive  class  of  ions  in 
various  gases  under  a  broad  range  of  conditions.  The  detailed  applica¬ 
tion  to  disturbed  atmospheres  requires  quantitative  data  on  competing 
ion-alectron  recombination  processes,  as  well  as  the  pertinent  rates  for 
attachment,  detachment,  and  ionizing  processes  under  the  atmospheric 
conditions  of  interest.  The  results  of  the  present  study  assume  added 
importance  in  the  study  of  disturbed  atmospheres  since  the  concentration 
of  polar  molecules  may  be  enhanced  significantly  by  chemical  and  ion- 
molecule  reactions  under  disturbed  conditions. 

In  studies  of  recombination  following  an  atmospheric  disturbance, 
it  has  become  clear  that  atomic  and  molecular  processes  which  thermalize 
the  original  supra-thermal  electrons  are  important,  especially  in  a 
nearly  neutral  atmosphere.  Elastic  electron-atom  and  electron-molecule 
processes  are  very  inefficient  for  the  critical  energy  region,  and  in¬ 
elastic  events  due  to  electronic  excitation  are  improbable.  Rotational 
excitation  of  polar  molecules  represents  an  efficient  energy  loss  mech¬ 
anism  resulting  in  the  thermalization  of  supra-thermal  electrons,  due  to 
the  very  long-range  electron-polar  molecule  interactor  (which  leads  to 
large  inelastic  scattering  cross  sections).  Since  the  effective  rates 
for  recombination  and  attachment  are  much  larger  for  thermal  as  opposed 
to  supra-thermal  electrons,  the  rate  of  thermalization  of  supra-thermal 
electrons  may  represent  the  controlling  processes  in  recombination  fol¬ 
lowing  a  severe  atmospheric  disturbance,  especially  since  the  disturbed 
atmosphere  may  have  an  enhanced  concentration  of  polar  molecules. 

Detailed  calculations  have  been  performed  to  derive  the  cross  sec¬ 
tion  for  rotational  excitation  of  CO  and  CN  by  electron  impact  in  Sec- 
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tion  III.  Angular  distributions  and  momentum  transfer  cross  sections 
have  been  obtained.  Since  the  elastic  scattering  cross  section  at 
these  energies  is  small  compared  to  the  inelastic  cross  section,  the 
full  implementation  of  a  close-coupling  formalism  is  required.  By 
studying  the  process  in  various  polar  molecules,  it  is  possible  to  de¬ 
termine  the  dependence  of  the  resultant  cross  sections  on  such  parameters 
as  the  dipole  moment  and  rotational  energy  separation  of  the  polar 
molecules.  Calculations  should  also  be  performed  for  water  vapor  and 
other  polar  molecules  of  atmospheric  interest;  however,  the  results  pres¬ 
ently  available  indicate  that  the  process  is  important  for  the  under¬ 
standing  of  the  behavior  ot  the  disturbed  atmosphere.  The  cross  sections 
are  also  of  importance  for  the  calculation  of  three-body  ion-electron 
recombination  in  polar  gases  as  discussed  above. 

In  studies  of  a  low  temperature,  high  neutral  density  atmosphere 
under  either  normal  or  disturbed  conditions,  significant  concentrations 
of  negative  ion  species  exist.  Ion-molecule  reactions  and  chemistry 
often  are  seen  to  lead  to  those  negative  ions  with  the  largest  electron 
affinity.  Investigations  of  the  loss  rates  of  negative  ions  are  required 
for  the  satisfactory  prediction  of  the  rate  of  decrease  of  atmospheric 
ionization.  Under  the  thermodynamic  conditions  usually  encountered, 
direct  collisional  detachment  has  a  very  small  rate  coefficient  because 
of  the  endothermic ity  requirements  imposed  by  the  relevant  electron 
affinities . 

Under  a  previous  program  (Ref.  1-3),  a  study  has  been  performed  of 
the  mutual  neutralization  process 


X+  +  Y" 


X  +  Y 


which  represents  one  of  the  most  important  mechanisms  for  the  destruc¬ 
tion  of  negative  ions  in  the  atmosphere.  Although  serious  discrepancies 
exist  between  the  theoretical  results  and  available  experimental  data, 
which  will  not  be  discussed  here,  all  studies  show  that  a  very  large  rate 
coefficient  exists  for  the  mutual  neutralization  process  at  thermal 
energies. 
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In  an  essentially  neutral  atmosphere,  such  as  that  encountered  in 
the  D-region,  the  associative  detachment  process 

X"  +  Y  -  XY  +  e 


may  be  as  important,  or  even  more  important  than  the  mutual  neutraliza¬ 
tion  processes  for  the  destruction  of  negative  ions  as  discussed  in  Sec¬ 
tion  IV  in  greater  detail.  For  atmospheric  species,  the  prediction  of 
the  pertinent  associative  detachment  cross  sections  constitutes  a  rela¬ 
tively  severe  theoretical  and  experimental  problem.  Quantitative  calcu¬ 
lations  have  been  performed  only  for  the  simple  prototype  system  H  +  H  . 
The  temperature  dependence  of  the  process  is  still  somewhat  uncertain, 
and  quantitative  predictions  require  very  accurate  potential  energy  curves 

and  coupling  matrix  elements,  ’n  order  to  study  the  process,  improved 

-  2 

accurate  potential  energy  curves  have  been  calculated  for  the  0~  (  JI  ) 

3  -  1  ® 

state  relative  to  the  02(X  curve  of  Schaefer-Karris  (Ref.  1-4). 

Certain  errors  detected  in  the  Schaefer -Harris  calculations  have  re¬ 
quired  a  recalculation  of  our  0 2  curves.  From  these  data,  accurate 

values  for  the  internuclear  separation  associated  with  the  crossing  of 
-  2  3  - 

the  C^C  Ilg)  and  02(X  curves  can  be  achieved  and  an  upper  limit  for 

the  rate  of  reaction  for  the  important  atmospheric  associative  detach¬ 


ment  process 


0  +  0  -  02  +  e 


can  be  obtained.  Although  the  accurate  ab-initio  calculation  of  poten¬ 
tial  energy  curves  for  associative  detachment  processes  in  atmospheric 
species  is  relatively  expensive  and  slow,  the  results  of  such  studies 
are  required  in  order  to  understand  the  process,  especially  since  the 
required  non-radiative  coupling  matrix  elements  can  seldom  be  determined 
by  other  methods.  The  necessary  data  on  the  required  potential  energy 
curves  cannot  usually  be  obtained  by  the  Rydberg-Klein-Rees  analysis  of 
experimental  spectroscopic  data. 

Under  the  current  program,  a  new  area  of  investigation  was  pursued 
involving  exploration  of  the  possibility  of  using  effective  potential 
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methods  to  calculate  molecular  wavefunctions  and  potential  energy  curves 
for  diatomic  molecules  at  intermediate  internuclear  separations  as  dis¬ 
cussed  in  Section  V.  Encouraging  results  were  obtained  for  the  simple 
system  chosen  for  the  initial  investigation,  However,  although 

these  preliminary  calculations  were  very  successful,  it  is  considered 
that  current  technical  limitations  inherent  in  the  methods  preclude  the 
derivation  of  accurate  data  at  the  present  time  for  most  systems  of 
atmospheric  interest  without  resort  to  an  extensive  and  moderately  ex¬ 
pensive  research  program. 

As  indicated  above,  the  four  processes  which  were  selected  for  de¬ 
tailed  investigation  under  the  present  program  are  discussed  individually 
in  the  following  sections  of  the  report.  Individual  contributing  inves¬ 
tigators  are  identified  in  each  section.  A  brief  conclusion  and  recom¬ 
mendations  for  future  study  discussion  is  presented  in  Section  VI. 
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II.  ELECTRON- ION  RECOMBINATION  IN  A  DENSE  MOLECULAR  GAS 
by  D.  R.  Bates 

ABSTRACT 

A  semi-quantal  method  is  dsveloped  for  treating  the  recombination 
of  electrons  and  positive  ions  in  a  dense  molecular  gas.  Extensive  cal¬ 
culations  are  carried  out  relating  to  recombination  in  hydrogen,  in 
nitrogen,  in  carbon  dioxide  and  in  damp  mixtures  of  gases.  This  latter 
case  is  of  special  interest,  recombination  being  greatly  facilitated  by 
the  presence  of  quite  a  small  amount  of  water  vapor.  A  very  high  re¬ 
combination  coefficient  is  possible.  For  example  in  a  nitrogen  +  hydro¬ 
gen  mixture  which  contains  207.  by  volume  of  water  vapor  and  which  is  at 

atmospheric  pressure  and  at  2000  K  the  recombination  coefficient  is 

-9  3  -1  -9  3  -1 

predicted  to  be  between  about  5.1  x  10  cm  s  and  8.7  x  10  cm  s  : 

here  the  lower  limit  is  calculated  by  the  semi-quantal  method  assuming 

that  the  energy  levels  of  the  recombining  system  are  hydrogenic  while 

the  upper  limit  is  calculated  by  a  purely  classical  method  assuming  that 

the  energy  levels  of  the  recombining  system  form  a  continuum. 

Recombination  coefficients  have  been  deduced  from  measurements  on 
the  ionization  in  flames  at  atmospheric  pressure.  The  values  obtained 
lie  satisfactorily  between  the  predicted  limits. 

A.  INTRODUCTION 

Three-body  ion-electron  recombination  processes 

X+  +  E  +  M-*X+M 

may  be  an  important  recombination  process  in  the  D  region  at  its  lower 
temperatures  and  higher  densities,  especially  for  those  ions  which  do 
not  undergo  chemical  reactions.  In  earlier  studies  it  was  shown  that 
much  higher  recombination  coefficients  were  obtained  if  the  third  body 
M  in  the  above  equation  exhibits  low  energy  degrees  of  internal  freedom 
as  is  the  case  for  a  large  class  of  molecular  third  bodies.  In  the 


Preceding  page  blank 


..^^35^y“0-~.  tv.!~  -,^!~v.> .  j-.TT.I,  tTw  —;■.  >  -V 


present  report,  attention  has  been  focused  on  the  polar  molecule  t^O 
as  the  third  body.  Calculations  have  been  carried  out  for  recombination 
in  several  reactions  of  diatomic  third  bodies  with  water  vapor.  The  re¬ 
sults  have  been  obtained  at  several  temperatures  and  mixing  ratios. 

For  an  ion  of  mass  23  a.m.u.,  a  recombination  coefficient  of  the  order 
-5  3-1 

of  10  "  cm  sec  can  be  obtained  for  ~  20 7,  by  volume  mixtures  of  1^0 
in  N2  at  250°K,  where  the  recombination  coefficient  is  expressed  as  an 
effective  two-body  coefficient. 

B.  BACKGROUND 

Recombination  of  electrons  and  ions  in  a  dense  neutral  gas  was  first 
treated  by  Thomson  (1924)  using  a  modification  of  his  model  for  three- 
body  ion- ion  recombination  (Massey  and  Burhop  1952).  An  entirely  dif¬ 
ferent  approach  was  introduced  by  Pitaevskii  (1962).  Taking  the  neutral 
gas  to  be  monatomic  he  showed  that  if  the  temperature  is  sufficiently 
low  the  recombination  process  may  be  regarded  as  classical  diffusion 
in  energy  space;  and  from  the  Fokker-Planck  equation  he  hence  derived  a 
simple  formula  for  the  recombination  coefficient.  Pitaevskii's  classical 
diffusion  method  was  later  applied  to  a  molecular  gas  by  Dalidchik  and 
Sayasov  (1966,  1967).  Its  main  defects  are  the  assumption  that  the 
negative  energy  levels  of  the  electron  in  the  field  of  the  ion  form  a 
continuum,  and  the  lack  of  explicit  allowance  for  radiative  transitions. 
Bates  and  Khare  (1965)  have  developed  a  semi-quantal  method  without 
these  defects.  However,  they  confined  their  attention  to  the  case  of 
a  monatomic  gas  where  all  collisions  are  elastic.  Their  results  agree 
with  those  of  Pitaevskii  only  at  very  low  temperatures.  It  is  of  interest 
to  extend  the  method  to  cover  recombination  in  a  molecular  gas.  Inelastic 
collisions,  involving  rotational  or  vibrational  transitions,  are  here 
important. 

The  process  with  which  we  shall  be  concerned  is  of  importance  in 
connection  with  ionization  in  flames. 

C.  THEORY 

We  shall  treat 

X++e+Z  -  X+Z  (II-l) 

10 


where  X+  is  an  atomic  or  molecular  ion  and  where  Z  is  a  neutral  molecule. 
In  this  section  we  shall  introduce  the  notation  and  quote  the  relevant 
formulae  of  the  semi-quantal  method. 

Averaging  over  the  initial  rotational  and  vibrational  levels  of  Z, 
which  we  suppose  to  be  in  thermodynamic  equilibrium  at  temperature  9, 
and  summing  over  all  possible  final  levels,  let  Jc(n,m)  be  the  rate  coef¬ 
ficient  for 

X(n)  +  Z(averaged)  -*  X(m)  +  Z(sunmed)  (11-2) 

n  and  m  being  the  principal  quantum  numbers  of  excited  (and  assumed  by 
hydrogenic)  initial  and  final  states  of  the  system  formed  by  recombina¬ 
tion.  Put 

e(»)  *  (11-3) 

where  aQ  js  the  atomic  unit  of  length.  Write 

<€{q)  -  A*(Z)  £  *  £  1(l  +  a+<H?+«)*'c(g+*.?-l-t)«xpW?+*)3.  (II-4) 

€-0  <-o 


where  N(Z)  is  the  number  density  indicated  and  q^  is  the  principal  quan¬ 
tum  number  of  the  highest  level  whose  population  may  rise  sufficiently 
to  deplete  the  population  of  the  continuum  appreciably  (which  number  may 
in  practice  be  taken  to  be  3  without  causing  appreciable  error);  write 


J*(?)  =  gM(?)exp[e(q)], 


Let  0^  be  the  radiative  recombination  coefficient.  From  the  work  of 
Bates  and  Khare  (1965)  we  have  that  the  rate  coefficient  describing  re¬ 
combination  due  to  the  combined  effect  of  collisions  and  radiative  transi¬ 
tions  is  closely 


<*cn  *  aCR+aR* 


(11-5) 


in  which 


v\  v  . 


(II-6) 


(II-7) 


The  atomic  mass  of  the  ion  and  the  chemical  symbol  for  a  molecule  of  the 
ambient  gas  will  be  shown  in  brackets  after  0^,R  where  this  is  judged 
desirable. 

Dissociative  recombination  also  occurs  if  the  ions  are  molecular- 
The  rate  coefficients  are  approximately  additive. 

The  transition  probabilities  required  for  the  radiative  conductances 
jg[q)  are  known  (see  Refs,  cited  by  Bates  and  Khare  (1965),  and  also  Wiese, 
Smith  and  Glennon  (1966)  and  Wiese,  Smith  and  Miles  (1969)).  Table  II*  1 
gives  the  values  of  j /(q)  for  atomic  hydrogen,  which  we  denote  by 
together  with  the  factors  f(n)  by  which  these  radiative  conductances  would 
be  reduced  if  transitions  to  level  n  or  lower  were  excluded.  The  informa¬ 
tion  contained  enables  estimates  for  complex  atoms,  which  are  accurate 
enough  for  present  purposes,  to  be  made  quickly. 
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TABLE  II- 1.  Radutiv*  cokdttctahom 
temperature,  6j  K 


ft 

250 

500  1000 

J*(?!H)/s-» 

_ _ _ A _ 

2000 

ft  01 

flW 

9 

» 

1.2* 

8.0“ 

1.7“ 

2.8“ 

0 

JW 

4 

8.9“ 

l.V» 

G.4U 

3.9“ 

0.51 

0 

6 

1.7** 

6.1  • 

1.0» 

4.4* 

0.88 

0.38 

6 

5.8“ 

8-1 1  ’ 

1.0“ 

1.1* 

0.73 

0.50 

7 

3.8“ 

5.7»- 

2.3* 

4.8* 

0.78 

0.58 

8 

1.3“ 

9.8* 

7.9* 

2.3* 

0-79 

0.83 

1  indkwe  ghre  the  powers  of  10  by  which  the  entries  should  be  multiplied. 


The  calculation  of  the  binary  rate  coefficients,  K(n,m),  and  hence  of 
the  collisional  conductances,  *(q)  ,  is  the  main  task. 


1 .  Rate  Coefficient  *(n,m) 

Consider  collisions  between  free  electrons  and  molecules  which  are 
in  a  specified  level.  Let  K^CE^jE^JdEj  be  the  rate  coefficient  for  col¬ 
lisions  in  which  the  energy  of  an  electron  is  changed  from  E,  to 
between  and  E2  +  dE^  and  the  internal  energy  of  a  molecule  is  in¬ 

creased  by  an  amount  I. 


Denote  the  masses  of  the  ion  and  electron  by  and  mg  respectively. 


Take  V  ,  and  V  _  to  be  the  velocities  of  the  electron  relative  to  the 
el  e2 


ion  before  and  after  the  collision;  and  take  and  V  to  be  the  cor¬ 


responding  velocities  of  the  molecule.  Write 


«  -  »a  _  -v  •» 

We  find  that  to  a  close  approximation 


(II-8) 


(H-9) 


where 


J  =  m,C7.(F„-Frt) 


(II-10) 


(Il-il) 


and 


S  “  (njlM-rr)  FM  •  (Km ~  Kn)* 


M^  being  the  reduced  mass  of  the  ion  and  molecule. 

In  the  first  instance  ve  shall  suppose  that  6,  which  is  in  general 
very  small  compared  with  I  +  J,  may  be  neglected.  We  shall  assume  that 


~-Vn «  U. 


(H-12) 


For  simplicity  we  shall  also  assume  that  the  e— Z  scattering  is  isotropic 
so  that  the  orientations  of  and  V '  ^  with  respect  to  U  are  random.  In 
many  cases  the  scattering  is  non-isotropic  (Bardsley  and  Reed  1968)  but 
the  error  arising  is  unlikely  to  be  serious. 


Represent  the  scattering  cross-section  and  the  normalized  classical 
distribution  functions  (Bates  and  Khare  1965)  by  a  (Ve|),  g(V£l)  and 
f(U)  respectively.  Put 


c.f 


0.  *>.#- 


(11-13) 


We  then  have  that 


KA«v  S,)d E,  -  iff Vncrt{V,l)g(V,1)f(U)dftlAiltdUdVn,  (11-14) 


where  S  signifies  that  the  integration  is  carried  out  over  all  accessible 
regions  in  which  the  final  kinetic  energy  of  the  electron  is  within  an 
interval  dE£  around  the  value  E2-  The  lowest  contributing  value  of  is 


F,  -  (27/m.)*  (/  >  0) 


>0)1 

0  (/  <  0)  j 


(11-15) 


and  the  highest  is  a  cut-off  value  V  defined  by 
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s 

1 
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4 
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fe 


P 


K 


K 


I 


K 


L 


(H-16) 


Tg  being  the  closest  the  electron  can  be  to  the  ion  and  yet  be  considered 


as  colliding  with  the  molecule.  For  a  given  value  of  V^,  it  is  clear 


that 


Ue  <  U  <  co, 


(11-171 


where 


Uc-  |  J !  /**Ki  +  ( -  2//«,)*]. 


(11-18) 


It  is  now  convenient  to  denote  the  greater  of  V  ,  and  V  -  by  V  >  and  the 

el  e2 

lesser  by  V<  and  to  put 


M>  •?>.<>,  m<»K^. 


(H-19) 


If  J,  U,  V>  and  V<  are  held  fixed  while  pi<  is  varied  the  extent  of  the 
range  u>  is 

at  -  2^,  0  <  |J|  <  me  U(V._  - F<) 


-  1 .  »>.  -  f<)  <  \j\  <  «,  UiV, +F<). 

meC/  r  > 


(11-20) 


We  find  that 


Kj(Ex,Et)  * 


•i 


(H-21) 


with 


™-<ZvXSU-'miL' 


(11-22) 


Set 


-  EJkO  *  A,,  -  EfjkB  m  A„ 
I]kO  =  >\  J  !k6  mi]. 
m*VUUdmr,  (««/*0ro)-A,  -  *c. 


(11-23) 

(11-24) 

(11-25) 
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and 


-  {*- v)»}*  9  X",  n'lMzk + (*-  v)i}»  9  x+, 
iMxzImt  9  l 

Kt{Ex,Et)dEt  s  ^(Aj,  Aj) dAj. 


By  straightforward  analysis  we  deduce  from  (11-21)  that 


^  =  L.  ^+^ijr(x)dx' 

X(x)  =  (*,  +  i,)  exp  [  -  X+]  -  (ij  -  it,)  exp  [  -  $iX~) 


where 


+  M  (g^)}  V[i«+)»-erf[i/X^, 


(11-26) 

(11-27) 

(11-28) 

(11-29) 


(11-30) 


(H-31) 

(11-32) 


I 

I 

i 


•S' 

I 

-i 


in  which 


X09J',  3  */(*-„)*, 


if  v  is  positive;  and  in  which 

xo  s  s  zi,  t,  ■  x*/(x-y)t, 


(11-33) 


(11-34) 


if  v  is  negative. 

As  a  partial  check  it  may  be  noted  that  if  we  take  xc  to  be  infinite 
we  may  deduce  from  (11-30)  that  the  rate  coefficient  describing  the  excita¬ 
tion  of  the  molecule  is 
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^^^^WfrTir  T, 


v*»  t-  -  v  ■•  ^-r»  ."j~"v ~^<t,^. <  7.^  ’„  ' 


fVwKiM^dF., 

J  v. 


(II-35) 


as  of  course  it  should  be. 


We  may  now  obtain  an  approximation  to  the  rate  coefficient  K(n,m|i,j) 
for  those  (X+  +  e)  -  Z  collisions  in  which  the  electron  changes  from  an 
orbit  or  principal  quantum  number  n  to  one  of  principal  quantum  number  m 
and  in  which  the  molecule  changes  from  level  i  of  energy  eJc©  to  level  j 


of  energy  e^k0  with 


e,-e,  =  v. 


(11-36) 


Following  the  procedure  adopted  by  Bates  and  Khare  (1965)  we  have  that 


(II-37) 


where  if  we  put 


At  s  e*/2«*a0W,  A,  =  e*/2m*ooi0, 


(11-38) 


(11-39) 


1 


K'Jn,mli,j)  =  exp  (A,— Ax — v)  x*(m,  »|  j,  t),  (11-40) 


and  being  the  statistical  weights  of  the  levels  indicated. 

Averaging  over  the  initial  levels  i  and  summing  over  the  final  lev¬ 
els  j  we  see  that  the  rate  coefficient  sought  is 


*{w,to)  «  K0(n,m)  +  Kx(n,m), 


(11-41) 


where  Kg(n,m)  is  '.he  contribution  from  elastic  scattering  (see  Bates  and 
Khare  1965)  and  where 

*The  simple  expression  which  Bates  and  Khare  gave  for  k.Q(n,m)  may  be 
recovered  from  our  express’'  lor  »>(n,m|i,j)  by  carrying  through  some 
elementary  analyses. 


(11-42) 


rx(n,m)  -  j-jr  2  S  *><exp(-e<)*(n,ro|*>j-/ 
JVP)  i  i+i 


is  the  contribution  from  inelastic  e-Z  scattering,  ^2C®)  being  the  internal 
partition  function  of  the  molecule.  We  shall  denote  the  contributions 
to  K^(n,m)  from  rotational  and  vibrational  transitions  by  K^(n,m)  and 
Ky(n,m)  respectively.  The  letters  affixed  as  subscripts  to  jc's  will  also 
be  affixed  to  other  quantities. 


Examination  of  equations  (11-30)  to  (11-32)  shows  that  the  rate  coef¬ 
ficients  are  sensitive  to  the  net  change  in  the  internal  energy  of  the 
complex  consisting  of  the  recombining  system  and  the  third  body.  Antici¬ 
pating  the  results  of  the  detailed  calculations  we  illustrate  this  in 
Table  11-2.  As  would  be  expected  physically  the  resonance  effect  exhibited 
becomes  sharper  as  the  temperature  6  is  decreased  or  as  the  reduced  mass 
“xz  is  increased.  Irregularities  in  the  variation  of  the  collisional 
conductance  *(q)  with  q  are  evident  in  some  of  the  tables  to  be  presented 
later.  These  are  due  to  resonances. 


Method  of  Dalidchik  and  Sayasov 


According  to  the  classical  diffusion  method  of  Pitaevskii  (1962) 
the  recombination  coefficient  is  given  by 


r  _  e» /tt\I  r  f°  |S|*exp(£/i0)dEl~l 
4  vH?/  [ J  J 


(11-43) 


where  &{E)  is  the  mean  of  the  square  of  the  change  in  the  energy  E  of  the 
electron  per  unit  time. 


Pitaevskii  considered  only  elastic  collisions.  Taking  the  scattering 
to  be  isotropic  and  the  cross-section  to  have  the  constant  value  a  he 
demonstrated  that 


(11-44) 
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T/  iLE  II-2.  Resonance  intoi  in  rotational  transitions 

OT  HYDROGEN 


rotation*! 


ion  maas/a.m.u. 


perature, 

6fK 

250 

quantum 

number, 

T 

2 

4  8 

*r(20, 1 1 J  J,  J + 2)/cm*  a-1 

_ A. _ 

co 

■ 

V 

0 

7.7--* 

4.5-“ 

3.1-“ 

1.9'“ 

1 

1.1-“ 

1.0-“ 

9.6-“ 

8.8'“ 

2 

2.0'“ 

1.3-“ 

9.6-“ 

6.3~“ 

.700 

0 

1.4-“ 

9.9-“ 

7.6-“ 

5.3"“ 

1 

1.1-“ 

1.1-“ 

1.1-“ 

1.0'“ 

2 

8.4-“ 

4.6-“ 

3.7-“ 

2.7-“ 

1000 

0 

2.5*“ 

1.9-“ 

1.5-“ 

1.2~“ 

1 

1.0-“ 

1.1-“ 

1.1-“ 

1.1-“ 

2 

1.4-“ 

1.1-“ 

9.3-“ 

7.6-“ 

2000 

0 

4.1-“ 

3.3-“ 

2.8  “ 

2.3-» 

1 

8.6-“ 

9.3~“ 

9.7-“ 

1.0-“ 

2 

2.6*“ 

2.1'“ 

!.8-“ 

1.6-* 

Note*.  The  internal  energy  lost  by  the  recombining  system  in  the  (20, 11)  transition  is  9101: 
where  k  in  the  value  of  the  Boltzmann  constant  (per  kelvin);  the  internal  energy  gained  by 
the  H,  molecule  in  the  (J,J  +  2)  transition  is  525k,  8751:  and  1225k  for  J  =  0,  1  and  2.  The 
indices  give  the  powers  of  10  by  which  >'<9  entries  should  be  multiplied. 
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Dalidchik  and  Sayasov  (1966)  showed  how  to  take  rotational  transitions 
into  account.  For  homopolar  molecules  they  used  the  Born  approximation 
to  the  cross-section  concerned  which  has  a  very  simple  form  (Gerjuoy  and 
Stein  1955).  Elementary  analyses  then  yields 


*(£)-#(  Z){ 


45(2m,)i  J’ 


(11-45) 


where  B  is  the  rotational  constant  of  the  molecule  and  geaj  is  the  quadru- 
pole  moment.  The  contribute ons  ^(E)  and  ^R(3)  to  &{ E )  are  of  course  addi¬ 
tive.  Dalidchik  and  Sayasov  (1967)  also  treated  linear  molecules  having 
a  permanent  dipole  moment.  Taking  the  cross-section  for  rotational  exci¬ 
tation  to  be  as  given  by  the  Born  approximation  (Takayanagi  1966)  they 
found  that 


(11-46) 


deog  being  the  dipole  moment.  Where  necessary  we  shall  distinguish  be¬ 
tween  ^R(JE)of  (11-45)  and  (11-46)  by  adding  the  letters  q  and  d  as  subscripts. 

Comparing  (11-44),  (11-45)  and  (11-46)  it  is  seen  that  collisions  involv¬ 
ing  rotational  transitions  are  predicted  to  become  increasingly  more 
effective  than  elastic  collisions  as  the  energy  of  the  electron  is  re¬ 
duced.  This  prediction  is  in  conflict  with  wliat  mignt  be  expected  from 
simple  physical  considerations  and  with  the  results  of  the  semi-quantal 
theory.  The  error  in  the  prediction,  which  is  not  of  great  importance  in 


*This  is  twice  the  value  given  by  Dalidchik  and  Saya< ov  (1966)  who 
omitted  to  allow  for  level  J  being  (2j+l)-fold  degenerate. 
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practice,  arises  because  of  the  treatment  of  Dalidchik  and  Sayasov  does 
not  take  fully  into  account  the  existence  of  a  threshold  energy  for  rota¬ 
tional  excitation. 

The  mean  of  the  square  of  the  change  in  E  per  unit  time  due  to  vibra¬ 
tional  transitions,  ,  was  not  discussed  by  Dalidchik  and  Sayasov. 

There,  is  unlikely  to  be  a  simple  analytical  expression  for  it  which  is 
accurate  enough  to  be  useful.  However,  if  the  relevant  excitation  cross 
sections  and  energies  are  known  of  course  be  obtained  by  carry¬ 

ing  through  the  necessary  integration  numerically. 

3.  Calculations 

Calculations  were  carried  out  relating  to  recombination  in  the 
presence  of  hydrogen,  nitrogen,  carbon  dioxide  and  water  vapor.  All  con¬ 
tributing  processes  in  hydrogen  and  nitrogen  were  investigated.  Atten¬ 
tion  was  however  concentrated  on  the  contribution  from  vibrational  transi¬ 
tions  in  the  case  of  carbon  dioxide  and  on  that  from  rotational  transitions 
in  the  case  of  water  vapor.  These  are  of  interest  in  connection  with  re¬ 
combination  in  gas  mixtures. 

We  shall  first  briefly  indicate  the  assumptions  made  about  the  rele¬ 
vant  molecular  properties. 

The  structural  parameters  were  taken  from  the  books  by  Herzberg  (1945, 
1950)  or  from  references  cited  in  them.  The  quadrupole  moment  of  is 
0.393eaQ  (Harrick  and  Ramsey  1952)  and  that  of  is  0.96ea^  (Smith  and 
Howard  1950);  the  dipole  moment  ot  ^0  is  0.736eag  (McClellan  1963). 

The  elastic  scattering  cross-sections  of  Hg  and  towards  slow  elec- 

o  2 

trcns  were  put  equal  to  jOa^  and  15Sq  respectivelj  (Moiseiwitsch  1962). 
Elastic  scattering  by  HjO  was  disregarded  since  scattering  by  a  rotating 
dipole  is  predominantly  inelastic  (Crawford,  Dalgarno  and  Hays  1967). 

The  values  adopted  for  the  rotational  excitation  cross-sections  of 
Hj  and  ^  were  those  given  by  the  distorted  wave  approximation  (Takay- 
anagi  and  Geitman  1965).  In  the  case  of  Hj  these  are  much  larger  than 
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the  cross-sections  given  fay  the  3orn  approximation  (Gerjuoy  and  Stein 
1955),  but  in  the  case  of  Nj  the  two  cross-section3  do  r.ot  differ 
greatly.  The  rotational  excitation  cross-sections  of  1^0  were  obtained 
by  substituting  the  line  strengths  calculated  by  Cross,  Hainer  and  King 
(1944)  info  the  formula  which  Crawford  (196?)  derived  for  asymmetric  top 
molecules  using  the  Bom  approximation,  ilalidchik  and  Sayasov  (1967) 
took  1^0  to  be  a  linear  dipole  molecule  having  a  moment  of  inertia  equal 
to  the  smallest  of  the  three  principal  moments.  Calculations  which  we 
have  carried  out  on  the  collisional  conductances  show  that  this  is  quite 
a  good  approximation. 


The  cross-sections  for  the  vibrational  excitation  of  P_  and  C0„  were 

* 

taken  from  measurements  by  Schulz  (1964)  and  by  Booess  and  Schulz  (1968) 
respectively;  those  for  the  vibrational  excitation  of  Nj  were  taken  from 
the  theoretical  investigation  by  Chen  (1964).  Reliable  information  on 
the  vibrational  excitation  of  HgO  is  not  available.  Takayanagi  (1966) 
has  suggested  that  in  the  case  of  a  linear  molecule  the  cross-section  Is 


approximately 


<r{v,J-+v+l,J+l) 


(11-47) 


in  which  M  is  the  reduced  mass  of  the  vibrating  atoms,  Re  is  the  equilibrium 
internuclear  distance,  I^_  is  the  vibrational  quanta,  the  other  symbols 
have  their  usual  significance  and  all  quantities  are  in  atomic  units. 
Tentatively  adopting  this  formula  for  each  of  the  vibrational  modes  we 
rather  arbitrarily  set  M  equal  to  3  My  where  is  the  mass  of  the 
proton  and  took  the  rotational  energy  levels  to  be  as  in  the  model  used 
by  Dalidchik  and  Sayasov  (1967)- 

The  distance  rc  introduced  in  (11-16)  is  ill  defined.  We  took  it  to 
be  0.1  nm.  Fortunately  the  results  are  insensitive  to  the  value  chosen. 

Little  need  be  said  about  the  numerical  procedure  which  is  straight¬ 
forward.  Care  is  needed  in  carrying  out  the  quadrature  in  (11-31)  because 

Following  a  private  communication  from  Dr.  Schulz,  the  threshold  energy 
of  1.0  eV  originally  reported  was  modified  to  0.53  eV  but  the  overall  pat¬ 
tern  of  the  variation  of  the  cross  section  was  retained. 
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of  the  sharp  peak  exhibited  by  the  integrand.  Severe  cancellation  occurs 
between  the  error  functions  in  (II-32)  so  that  they  must  be  represented 
very  accurately.  We  used  the  22-term  even  Chebyshev  series  for  4x  erfx 
through  the  region  in  which  the  independent  variable  x  is  less  than  4.0 
and  used  the  9-term  even  Chebyshev  series  for  trfceexp{**)erfc*  elsewhere 
(Clenshaw  1963). 

Transitions  involving  levels  for  which  q  is  very  high,  say  above  sons 
value  ,  are  caused  mainly  by  elastic  collisions  (except  in  the  case  of 
polar  molecules).  It  is  sometimes  convenient  to  treat  them  by  the  method 
of  Pitaevskii  (1962).  Hiis  simplifies  the  evaluation  of  the  infinite 
summation  appearing  in  <XI-6):  thus  we  obtain 


V  M?) _ (  &  /*(?)  },p(a) 

-£«*{?)+.*(?)  “  W-^4i  «•(?)+ ^(9)1  ^ 


(11-48) 


where 


Jexpl  ag»o#w)}-  (II’49) 


Results  and  Discussion 


a .  Hydrogen 

Figures  n-1  and  IT-2  show  the  results  obtained  on  the  rec ona> inat ion  of 
an  ion  of  atomic  mass  2  in  molecular  hydrogen.  The  values  of  the  recom¬ 
bination  coefficient,  a^B(2,H,)  are  very  high.  They  greatly  exceed  those 
for  o^R(4,He)  as  calculated  by  Bates  and  Khare  (1965).  Thus  if  the  ambient 

IQ  *3  _*7  O  ■» 

number  density  is  3  x  10  cm  ,  “crI^.H,)  is  8  x  10  cm  s  at  250  K 

-9  3  -1 

and  is  8  x  10  '  cm  s  at  2000  K,  whereas  the  corresponding  values  for 

“8  3  -I  *11  3  *1 

a^(4,Hft)  are  5  x  10  cm  s  and  2  x  10  cm  s  .  There  are  two 
main  reasons  for  the  difference:  the  reduced  mass  of  the  ion  and  third 
body  is  less  in  the  first  case  than  in  the  socond;  and  the  cross-section 
of  an  H2  molecule  for  the  elastic  scattering  of  slow  electrons  is  greater 
than  that  of  an  He  atom.  As  the  temperature  is  raised  the  contribution 
to  0^(2, Ht)  from  rotational  and  vibrational  tiansitions  increases  in 
importance  (Figure  II-2). 


number  density,  Ar(H,)/cna 

Figure  1 1*  1  •  Recombination  coefficient  for  ZT + e+H,  -* 

molecule  number  density.  The  temperature  0(K)  is  ir 
of  the  ion  is  taken  to  be  2  n.m.u. 


temperature,  6[K 

Figure  II -2 «  Ttcoombination  coefficient  aj*  for  Z*  +  o+H,  -*•  Z+H,  as  function  of  the  tempera¬ 
ture.  The  hydrogen  molooule  number  density  A’(H,)  (cm-’)  is  mdioatod  on  each  set  of 
curve*.  The  mass  of  tlio  ion  is  taken  to  bo  2  a.tn.u.  Curve  o,  all  types  of  transition 
included ;  curve  6,  vibrational  transitions  oxoludcd;  ourvo  c,  vibrations]  and  rotational 
transitions  excluded. 


Table  11-3  provides  information  on  the  recombination  mechanism  in 
the  form  suggested  by  Bates  &  Khare.  It  gives  (i)  the  critical  principal 
quantum  number  q  at  which  the  term 

A 


t  „  Mia) 

*  Vi9)  +  *(q) 


(11-50) 


of  the  summation  in  (n-7)  is  greatest;  (ii)  the  values,  ?ix-  and  fyx*of  q 
below  and  above  at  which  t(q)  is  closest  to $<(9x) ;  (iii)  the  actual 
ratios 


r*x-  =  HHx-Mlx).  rix.  «  t(9ix+)/<(5x). 


(11-51) 


and  (iv)  the  ratios  “  1)»  P(?x*“l)  and  Mfyx*  —  l)where  4(a)  is  as  defined  in 

(II-6)  sc  that  the  extent  to  which  these  ratios  fall  below  unity  indi¬ 
cates  the  importance  of  radiative  transitions.  Rotational  and  vibrational 
transitions  make  q^  larger  than  it  is  for  helium.  Moreover,  they  keep 
radiative  transitions  insignificant,  whereas  in  the  case  of  helium  such 
transitions  have  a  marked  effect  at  the  lower  densities  and  high  tempera¬ 
tures  studied. 

For  comparison  some  calculations  were  carried  out  by  the  classical 
diffusion  method.  A  trifling  modification  was  made  to  the  formula  of 
Dalidchik  and  Sayasov  (1966) :  the  evaluation  of  &r ,(E)  was  carried  through 
without  the  simplifying  assumption,  which  leads  to  (11-45),  that  the  rota¬ 
tional  constant  B  is  very  much  smaller  than  1:8,  this  assumption  being 
invalid  at  the  lower  temperatures  of  interest  in  the  case  of  hydrogen. 

In  fact  the  modification  had  little  influence  on  the  final  results.  The 
derived  values  of  some  of  the  recombination  coefficients  are  given  in 
Table  II-4  together  with  the  corresponding  values  obtained  by  the  semi- 
quantal  method.  To  avoid  complicating  the  comparison  these  last  are  also 
based  on  the  Born  approximation  to  the  rotacional  cross  sections  and  do 
not  include  the  contribution  from  vibrational  transitions.  As  would  be 
expected  the  classical  diffusion  method  gives  larger  recombination  coef¬ 
ficients  because  of  the  assumption  that  the  negative  energy  levels  of 
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TABLE  II-3  .  RECOMBINATION  MECHANISK  TOR  1098  OT  ATOXIC  MASS  2 

IK  MOLECULAR  HYDROGEN 


temperature, 

0/K 

10 

25 

>36 

125 

0.50 

~0.5 

1.0 

1.0 

1.0 

12 

(10,21) 

>30 

250 

0.50 

~0.5 

1.0 

1.0 

1.0 

8 

11 

21 

600 

0.66 

0.51 

0.07 

1.0 

1.0 

6 

8 

16 

1000 

0.43 

0.45 

0.05 

0.09 

1.0 

4 

7 

11 

2000 

0.58 

0.63 

0.94 

0.00 

1.0 

The  entries  for  each  temperature  are  almost  independent  of  the  number  density  N(Ht) 
provided  this  it  not  below  about  10”  cm-*.  They  ere  arranged  in  the  form 

<7ix-  Qix  *■ 


9ix- 

r\x- 

Mtlx—U 

ihe  Dotation  being  as  defined  in  the  text. 


?4x* 

ri*' 


TABLE  II -4.  Comparison  of  values  of  recombination  coefficient  obtained  by 
CLASSICAL  DIFFUSION  AND  SEMI-QUANTAL  METHOD  FOR  IONS  OF  ATOMIC  MASS 
2  nr  MOLECULAR  HYDROGEN 

temperature,  &IK 


125 

500 

2000 

10” 

(i) 

2.6-* 

5.5*” 

1.3*“ 

(») 

i.o-» 

2.8-” 

5.6*” 

10” 

(i) 

2.6-* 

5.5“ 

1.3“ 

(”) 

1.6“ 

2.0“ 

1.0*” 

10“ 

(») 

2.6-* 

5.5“ 

1.3“ 

(“) 

1.0“ 

1.5“ 

1.8-” 

Nolei.  (i)  indicates  cho  classical  diffusion  method  and  <u)  the  semi-quanta]  method  using 
seme  cross-section*  (tec  text).  Vibrntior.nl  transitions  arc  ignored. 

The  inmost  give  the  powers  of  10  by  winch  the  eutm«  should  bo  multipliod. 


an  ion  form  a  continuum.  The  effect  is  especially  marked  when  the 
density  and  temperature  are  high. 

The  results  in  Table  II-5  show  that  an  increase  in  the  mass  of 
the  ion  from  2,  on  the  chemical  scale,  to  infinity  causes  the  recombina¬ 
tion  coefficient  to  fall  to  about  30  to  40%.  Their  detailed  pattern  is 
rather  complicated  because  of  the  variation  with  temperature  of  the  rela¬ 
tive  importance  of  the  three  types  of  collision  involved.  A  significant 
increase  in  is  normally  accompanied  by  a  change  in  the  structure  of 
the  system  and  therefore  by  changes  in  the  radiative  and  collisional 
conductances  which  are  not  taken  into  account  in  Table  II-5.  The  radiative 
conductances  are  relatively  unimportant  (see  Table  II-3).  Owing  to  resonance 
effects  some  of  the  rate  coefficients  which  determine  the  collisional 
conductances  may  be  greatly  affected  by  the  differences  between  the  ener¬ 
gies  of  the  levels  of  X  and  the  assumed  hydrogenic  values.  The  recombina¬ 
tion  coefficient  is  however,  likely  to  be  affected  to  only  a  minor  extent 
since  it  depends  (see  Table  II-3)  on  a  considerable  number  of  transitions. 

An  upper  limit  to  the  recombination  coefficient  at  a  given  temperature  is 
provided  by  the  classical  diffusion  method  (provided  the  gas  density  is 
high  enough  for  radiative  transitions  to  be  unimportant).  The  last  row  of 
Table  II-5  shows  the  upper  limits  obtained  using  the  elastic,  rotational 
vibrational  cross  sections  adapted  in  the  calculations  by  the  semi-quantal 
method,  and  evaluating  SF-^E)  and  &v(E) by  numerical  quadrature.  The  factor 
by  which  the  upper  limit  exceeds  the  corresponding  recombination  coef- 


TABLE  11-5.  Influence  of  mass  Mx  of  ion  on  recombination  coefficient 

IN  MOLECULAR  HYDROGEN  :  N{ Hf)  =  101*  cm-* 


U)in)H'raturo,  0.  K 


-Vs  •'«  jn.u. 

2 

250 

500 

■*CR 

1000 

cm***1 

2000 

2.8-’ 

5.T-* 

1.2-* 

2.S-* 

4 

2.2”’ 

5.W* 

1.1 -» 

2.5~* 

8 

1.0-' 

4.5-* 

1.0-* 

2.4-* 

■0 

!.?-» 

4.2-* 

9.2-* 

2.0-* 

upper  limit  (cIamicaI 
diffiuion  method): 

.Vx  infinite  1.0-* 

2.7-’ 

4.0-* 

5.3- • 

Th®  indice*  giv®  the  power*  of  10  by  which  the  on  trio*  ithould  bo  multiplied. 
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ficient  decreases  as  the  temperature  is  increased.  This  is  because  the 
inhibiting  effect  of  lack  of  balance  between  the  initial  and  final 
internal  energies  of  the  colliding  systems  is  less  marked  at  high  tempera¬ 
tures  than  at  low  temperatures. 

The  recombination  coefficient  in  a  mixture  of  neutral  gases  is 
not  the  sum  of  the  recombination  coefficients  in  each  of  the  constituent 
gases.  However,  the  collisional  conductances  are  additive;  and  using 
formula  (II-7)  the  recombination  coefficient  may  readily  by  obtained  from 
these  and  the  radiative  conductances  (see  §  2  in  particular  Table  II-l). 

The  important  collisional  conductances  for  molecular  hydrogen 
are  given  in  Table  II-6.  Except  hear  a  resonance  they  are  decreasing  func¬ 
tions  of  the  i.  .1  mass. 


TABLE  II- 6.  Collisional  conductances  ix  HYDROOExt 


temperature,  O' K 


quantum 

level 

ft 

2i>J 

500  1000 

Vlq)lX{  H,)a-« 

4. _ 

2000 

V 

ion  moss  2a.m.u. 

X 

4 

4.8-* 

2.6-> 

8.3-* 

5.0'7 

5 

1.1*= 

1.8-= 

8.6** 

7.3"7 

6 

3.2*= 

4.3-* 

2.0-* 

0.4-7 

f 

3.0-= 

S.4-* 

9.8-7 

5.4 -7 

8 

1.9'* 

3.3-* 

9.’.-7 

6.0-7 

9 

1.2-* 

3.3-* 

i.o-* 

7.3  7 

10 

6.6'» 

3.0-* 

1.2* 

8.S-7 

11 

1.2-= 

2.2-* 

1.2-* 

1.0-* 

12 

1.1-* 

2.4-* 

1.4-* 

1.2-* 

13 

i:2-* 

2.5-* 

1.6-* 

1.5-* 

14 

7.2** 

2.C-* 

1.8-* 

i.s-» 

15 

7.0 -* 

2.8~* 

2.0 -* 

2.2-« 

16 

3.7“* 

2.H-* 

2.3-* 

2.6-* 

17 

6.9-‘ 

3.0-* 

2.7-* 

3.4-* 

18 

5.9-‘ 

3.3-* 

3.!  -* 

4.2-* 

19 

5.9-* 

3.6-* 

3.7* 

5.1-* 

20 

5.6-* 

3.9-* 

4.3-* 

0.2 -* 

21 

5.6-* 

4.3-* 

3.0-* 

- 

22 

5.8-* 

4.9-* 

— 

— 

23 

6.1-* 

— 

— 

— 

temperature,  6'K. 


quantum 

level 

250 

500  1000 

*"(?)/-Nr(H1)/s-' 

2000 

9 

4 

ion  mass  4  a.m.u. 

2.1-*  1.4-=  6.5-* 

4.6  7 

5 

5.1*' 

1.3-* 

7.0-* 

6.5  7 

6 

2.3-= 

3.7-* 

1.8-* 

6.2  7 

7 

2.5-’ 

7.3* 

9.2-7 

5.3- 7 

8 

1.G-* 

3.0-* 

9.0-7 

5.9  7 

9 

1.2-= 

3.3-* 

i.o-* 

7.1  7 

10 

5.6 -» 

3.0-* 

1.1-* 

8.3-7 

4 

ion  mass  8  a.m.u. 

3.9- «  9.4-*  5.5-* 

4.1  7 

5 

2.7-' 

9.2* 

5.9-* 

0.0  7 

6 

2.0-* 

2.7-* 

1.7-* 

c.o-7 

7 

0.7-* 

9.0-7 

5.3-7 

8 

1.4-* 

2.7-* 

S.9-7 

5.9 -7 

9 

1.3-* 

3.3-* 

i.o-* 

7.0-7 

10 

0.0-* 

3.0-* 

1.1-* 

8.0 -7 

4 

ion  mass  infinite 

1.1-'  5.C-*  4.3'* 

3.5  7 

•7 

1.0*' 

5.C-* 

4.6-* 

5.0  7 

(» 

1.8-= 

3.3-» 

1.6-* 

5.9  7 

7 

1.9-* 

0.4“ 

8.8‘7 

.">.3  7 

$ 

1.2-* 

O  «(-• 

9.0-7 

o.o 7 

9 

1.4-* 

4.0-* 

i.o-* 

7.0-7 

10 

0.3  » 

3.7-* 

1.1-* 

7.9 7 

t  More  extensive  tables  ore  pi  von  bv  Malaviya  (1970). 

Hie  indices  give  the  powers  of  10  by  which  the  entrii-s  should  be  multiplied. 
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Nitrogen  and  Carbon  Dioxide 

Results  of  calculations  on  Ocr(28|Nj)  by  the  semi-quantal  method, 
with  the  radiative  conductance  taken  to  be  f(2)s/(q\S)  of  Table  II-l  are  dis¬ 
played  in  Figure  11-3  along  lines  of  constant  temperature.  Figure  II-4  shows 
some  of  the  same  results  along  lines  of  constant  density  and  also  shows 
the  effect  of  neglecting  vibrational  transitions.  For  comparison  pur¬ 
poses  results  of  calculations  using  formulae  (11-43)  and  (11-45)  of  the  clas¬ 
sical  diffusion  method  are  included  in  this  figure. 

As  may  be  seen  collisional  radiative  recombination  in  nitrogen 
is  much  slower  than  in  hydrogen.  There  are  two  main  reasons  for  the  dif¬ 
ference.  First,  elastic  collisions  are  unimportant  in  the  case  of  nitro¬ 
gen  because  of  the  molecules  being  so  massive,  whereas  they  are  quite 
significant  in  the  case  of  hydrogen.  Secondly  the  rotational  constant  B 
is  considerably  smaller  for  nitrogen  than  for  hydrogen  which  renders  the 
rotational  transitions  less  effective  (see  Figures  II-2  and  11-4). 


JO'-"  ID"  10J 
number  density.  A'(K3)/cm"* 

Figure  Recombination  coefficient  «(.»  for  Z*  +  e+N,  -*Z  +  X,  as  function  of  the  nitrogen 

molecule  number  density.  The  temperature  6{K)  is  indicated  on  each  curve.  The  mass 
of  the  ion  is  taken  to  be  28  a.m.u. 
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temperature,  OfK 

Figure  II-4.  Rcoombination  coefficient  aj*  for  Z'  +  e+iv,  -vZ+N,«*  function  of  the  tempera¬ 
ture.  The  nitrogen  moleculo  number  density  .V(Nt)  is  indicated  on  eaoh  curve.  The  mass 
of  the  ion  is  taken  to  be  28  s.m.u.  Curve  a.  all  types  of  transition  inoiuded;  curve  b, 
vibrational  transitions  excluded  (full  line,  semi-quanta! ;  broken  line,  olsssioal  diffusion). 


The  results  obtained  by  the  classical  diffusion  method  differ 
markedly  from  the  results  obtained  by  the  semi-quantal  method  (Figure  II-4). 
This  is  because  the  main  resistance  to  the  down-flow  of  the  recombining 
electrons  arises  from  the  wide  gaps  (ignored  in  the  classical  diffusion 
method)  between  the  levels  of  low  principal  quantum  number.  These  gaps 
are  mainly  bridged  by  collisions  involving  vibrational  transitions  (Fig¬ 
ure  11*4) . 

Table  II-7  shows  how  the  recombination  coefficient  derived  using 
the  semi-quantal  method  depends  on  the  mass  of  the  ion.  If  is 
expressed  in  units  of  the  mass  of  the  molecule  the  dependence  is  rather 
similar  to  that  for  hydrogen  (Table  H-5).  To  regard  as  a  quantity  which 
may  be  varied  in  isolation  is  of  course  artificial.  The  remarks  made  in 
this  connection  towards  the  end  of  3  are  relevant.  Little  guidance 
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TABLE  11-7 .  Ijtrujxuci  or  mam  Mx  or  ion  on  mcokbination  conrncixNT 
IN  SIOLECULAB  NITROOXN :  N{Ht)  -  10“  cm*5 

temperature,  OfK 


Af_ /a  m  it 

250 

600  1000 
a^j cm*  *~l 

_ A _ 

2000 

7 

1.7** 

5.3"* 

1.8"* 

8.1"«* 

14 

1. 1~* 

3.7*» 

1.4"* 

6.6"1* 

28 

7.6* 

2.7*» 

1.1-* 

5.5"“ 

56 

5.8~l 

2.1-* 

9.4*“ 

4.7"“ 

112 

4.8'* 

1.8-* 

8.2"l# 

4.2-“ 

upper  limit  (classical 
diffusion  method): 

Mx  infinite  2.4'* 

3.8"7 

5.7"* 

7.2-» 

The  indices  give  the  powon  of  10  by  which  the  entries  should  be  multiplied. 


is  provided  by  the  upper  limits  obtained  from  the  classical  diffusion 
method.  These  are  given  in  the  last  row  of  Table  11*7.  Each  is  very  much 
greater  than  the  corresponding  value  of  the  recombination  coefficient 
calculated  using  the  semi-quantal  method.  This  is  mainly  because  of  the 
closeness  of  the  molecular  energy  levels  and  the  largeness  of  the  reduced 
mass . 

Collisional-radiative  recombination  in  pure  nitrogen  gas  is  too 
slow  to  be  of  truch  interest.  Instead  of  presenting  further  results  on  the 
rate  of  this  process  we  give  in  Table  II-8  some  collisional  conductances 
which  may  be  required  in  calculations  on  recombination  in  gas  mixtures 
.  ontaining  nitrogen. 

In  the  case  of  carbon  dioxide  collisional  conductances  have  been 
computed  for  the  lower  levels  where  vibrational  transitions  dominate. 

The  values  obtained  have  been  tabulated  by  Malaviya  (1970).  They  are 
smaller  than  the  corresponding  values  for  nitrogen  and  will  not  be  given 
here. 

c .  Water  Vapor  and  Damp  Gases 

Table  II-9  contains  a  selection  of  the  calculated  rotational  con¬ 
tributions  to  the  collisional  conductances  of  water  vaoor.  Most  are 


TABLE  II -6.  COLLIHONAL  CONDUCTANCES  IN  NITBOaENf 

temperature,  OfK  temperature,  0/ K 


quantum  260 
level 

9  , - 


600  1000 

WW,!/*-' 


ion  sum  7  tan. 


quantum  250 
level 

3  < — 


50C  1000 


ion  man  26  a.m.u.  (coni.) 


4 

2.4+l 

2.7-* 

2.3-* 

1.6-’ 

13 

6.4-’ 

2.1-’ 

1.6-’ 

2.0-i 

5 

8.9** 

4.4~* 

3.6“* 

2.7~* 

14 

5.6-’ 

2.1-’ 

2.0*1 

2.8-’ 

0 

2.4*1 

2.8-* 

2.7-* 

2.2-* 

10 

4.9-’ 

2.7-i 

2.9-’ 

4.4-’ 

7 

i.o-» 

1.3-* 

4.3-* 

7.4-1 

16 

9.1-i 

4.2-’ 

4.4-i 

6.6-1 

8 

2.2-* 

3.4-» 

1.S-’ 

I.6-1 

17 

9.0-’ 

5.7-’ 

6.4-’ 

9.2-i 

B 

2.2** 

3.4-i 

I.6-1 

l.O-i 

18 

1.0-* 

7.4-’ 

9.0-’ 

— 

10 

2,0-4 

3.7-i 

1.9-i 

I.8-1 

19 

1.2-* 

9.5-’ 

— 

— 

20 

1.3“* 

. 

_ , 

_ 

ion  man  I4a.m.u. 

4 

3.1-» 

3.6-* 

3.0-* 

2.7-* 

ion  maa*  66  a.m.u. 

6 

3.5-1 

4.0-* 

3.9-* 

3.9-* 

4 

4.1*1 

4.8“» 

4.2-* 

2.8-* 

6 

3.1-* 

3.6-* 

3.6-* 

3.6-* 

5 

4.1  -* 

5.0-* 

4.9'* 

3.8-* 

7 

1.2-» 

1.6-* 

5.1-* 

9.3-’ 

6 

4.1-* 

5.0-* 

4.8-* 

3.7-* 

8 

1.4-* 

2.8-’ 

1.4-i 

1.3-1 

4 

1.3-* 

1.9-* 

6.3-* 

1.0-* 

9 

1.4-* 

2.8-’ 

1.6-1 

I.6-1 

8 

5.5-1 

1.7-1 

l.l-» 

1.2-’ 

10 

1.4~* 

3.0-’ 

I.8-1 

1.7-i 

9 

5.5-1 

1.7-i 

1.2-i 

1.3-i 

ion  maia  28  a.m.u. 

10 

5.0-1 

1.9-» 

1.6-’ 

1.3-’ 

4 

3.6*1 

4.2-* 

3.0-* 

3.2-* 

6 

3.7i-i 

4.4-» 

4.3-* 

4.4"* 

ion  man  112  a.m.u. 

8 

3.7-i 

4.4-* 

4.2-* 

4.2-* 

4 

4.4*1 

6.2-* 

4.6-* 

3  1-* 

7 

1.3-i 

1.S-* 

5.8-* 

1.0-* 

6 

4.4-* 

6.4** 

5.2- 

4.1-* 

8 

3.6-1 

2.2-’ 

1.3“’ 

1.3-’ 

6 

4.4-» 

6.4-* 

6.2-* 

4.0-* 

9 

8.6"  ’ 

2.2-i 

1.3-i 

1.5-1 

7 

1.3-1 

2.0-* 

6.6* 

1.1-* 

10 

8.0-5 

2.4-1 

I.6-1 

1.0-’ 

8 

3.9-- 

1.4-’ 

1.1-’ 

1.1-’ 

11 

1.1-* 

3.3-1 

2.0-i 

1.9-i 

9 

3.9-i 

1.4-i 

1.1-’ 

1.3-i 

12 

1.3* 

3.6-’ 

2.1-’ 

2.2-i 

10 

3.9-i 

1.6-i 

1.4-’ 

1.5-’ 

t  More  extensive  table*  are  given  by  Malaviya  (1970). 

Ths  indie**  give  the  power*  of  10  by  which  the  entries  should  be  multiplied. 


several  powers  of  tea  greater  than  the  corresponding  total  collisional 
conductances  of  non-pol^  gases. 


Using  formula  (11-47)  the  vibrational  contributions  at  the  lower 
principal  quantum  numbers  were  computed.  They  were  found  to  be  only 
moderate  in  magnitude.  In  treating  gas  mixtures  like  damp  hydrogen  or 
damp  nitrogen,  it  is  a  sufficient  approximation  to  replace  the  vibra¬ 
tional  contribution  from  the  water  vapor  by  the  vibrational  contribution 
from  an  equal  quantity  of  the  major  gas. 
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TABLE  11*9.  COLUSIOXAL  CONDUCT AUC*a  IN  WATXB  VAPOU*t 


temperature,  6[K  temperature,  6/ K 


quantum  2  SO 
level 

500  1000 

V(9)/tf(HtO)/.-* 

A 

2000 

quantum  250 
level 

500  1000 

«’(g)/2^(H,0)/.-* 

2000 

¥ 

11  t 

2011  man  9 1  tq.u. 

ion  maae  36a.m.u. 

a 

8.6-1* 

6.2-*' 

4.1-** 

6 

3.4-** 

2.4-** 

3.5-** 

9.1->* 

7 

2.5-» 

2.9-* 

2.4-* 

6.0-* 

7 

2.4-* 

2.8-’ 

2.4-» 

6.1-* 

8 

1.6-* 

1.2-* 

8.0-* 

6.2-* 

8 

1.6-* 

9.5-* 

6.6-* 

5.8-* 

9 

1.6-* 

2.1-* 

6.6-* 

2.3-* 

9 

1.5-* 

2.1-* 

6.7-* 

2.6-» 

10 

4.8-* 

3.2-* 

7.2”* 

3.5-* 

10 

4.7-* 

3.3-* 

7.6-* 

3.8-* 

11 

2.0'* 

3.6-* 

i.c-* 

1.0-* 

11 

1.6-* 

3.2-« 

1.6-* 

i.o-* 

12 

3.4'* 

e.3-« 

2.7-* 

1.8-* 

12 

2.8-* 

6.8-* 

2.7-* 

1.8-* 

13 

2.9-* 

7.6“« 

4.0-* 

3.0 -* 

13 

2.7-» 

7.2-* 

4.1-* 

3.1-* 

14 

4.5-» 

1.1-* 

6.8-* 

6.5-* 

14 

4.6~* 

1.1-* 

6.8-* 

5.6-« 

16 

3  6-* 

1.4-* 

1.0-* 

9.0-* 

15 

3.5-> 

1.4-* 

1.0-* 

9.1-* 

ion  mean  18a.m.u. 

6 

3.3-** 

1.2-** 

1.1-** 

2.5->* 

7 

2.6-* 

2.9-’ 

2.4-* 

6.3-* 

8 

1.6-* 

i.o-» 

7.2“* 

6.0-* 

9 

1.6-* 

2.1-* 

6.7-* 

2.4-* 

10 

4.8-* 

3.3"* 

7.6-* 

3.7-* 

11 

1.7-* 

3.4-* 

1.6-* 

1.0-* 

12 

3.1-* 

6.1-* 

2.7-* 

1.S-* 

13 

2.8-* 

7.3"* 

4.1-* 

3.1-* 

14 

*  6-* 

1.1-* 

6.8-* 

5.6-* 

15 

3.5-* 

.4-* 

1.0-* 

9.1-* 

t  Morn  extensive  tables  are  given  by  Mala  viva  (1970). 

The  indioes  give  the  powers  of  10  by  which  the  entries  should  be  multiplied. 


Figure  II-5  shows  the  calculated  collisional  radiative  recombination 
coefficients  in  some  damp  gases  as  functions  of  the  percentage  by  volume 
of  water  vapor,  a  pressure  of  one  atmosphere  being  maintained.  As  may  be 
seen  the  presence  of  a  small  amount  of  water  vapor  increases  the  recombina¬ 
tion  coefficient  markedly.  The  increase  is  due  to  the  rise  in  the  col¬ 
lisional  conductances  through  levels  of  high  principal  quantum  numbers 
(see  Tables  11-3  and  11-10). 

The  classical  diffusion  method  predicts  faster  recombination 
(especially  at  the  lower  temperatures)  since  substituting  from  (11-46) 
into  (  11-43)  gives 


1.2  x 


10“**tV(HjO) 


cra^-*. 


(11-52) 


33 


iarnal 


^gure  II- 5 .  Recombination  coefficient  ate,  in  damp  gas  mixtures  at  atmospherio  pressure  as 
ranotion  of  the  percentage  by  volume  of  water  vapour.  The  temperature  8(K)  is  indicated 
on  each  set  of  curves.  The  mass  of  the  ion  is  taken  to  be  23  a.m.u.  Curve  1,  hvdrogen- 
0UT7O  2,  nitrogen;  curve  3,  hydrogen-nitrogen  -.2). 


% 


TABLE  II- 10.  RiOOMRIHATIOR  MKJHANISM  FOB  loss  OF  ATOMIC  MASS  23  IK  A  DAMP 
hitbogbh+hydbogeh  (n(Nt):n(Ht)  ::2: 1)  mixture  at  atmospheric  pressure 

water  vapour  (%) 


temperature, 

* 

- - A. 

32 

6  IK 

7 

8 

9 

1 

8 

9 

500 

0.0 

0.4 

0.1 

0.1 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

7 

8 

10 

0 

8 

9 

1000 

0.1 

0.5 

0.1 

0.1 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

t 

S 

9 

G 

7 

9 

2000 

0.4 

0.6 

0.3 

0.3 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

The  basic  cause  of  the  difference  may  be  seen  by  considering 
the  integrand  of  expression  (11-43)  for  the  recombination  coefficient. 

It  goes  through  a  maximum,  corresponding  to  the  conductivity  going  through 
a  minimum,  where  | E j  is  3k 0  which  is  where  the  principal  quantum  number 
q^  is  2300-1/2  if  the  energy  levels  of  X  are  hydrogenic.  However,  the 
semi-quantal  calculations  show  that  the  rotational  contribution  to  the 
collisional  conductance  does  not  begin  to  increase  as  q  passes  below 
2300-1/2;  instead  it  continues  to  fall  more  and  more  steeply  because  of 
the  widening  separation  between  the  energy  levels.  The  collisional  con¬ 
ductance  does  indeed  increase;  but  this  is  due  to  the  growth  in  the 
vibrational  contribution.  Formula  (11-52)  does  not  take  the  vibrational 
contribution  into  account  which  introduces  an  error  in  the  opposite  sense 
to  that  arising  from  the  overestimation  of  the  rotational  contribution. 
Since  the  two  errors  are  entirely  unrelated  they  can  scarcely,  in  general, 
approximately  cancel. 

The  rate  ciefficients  for  the  rotational  transitions  of  water 
greatly  exceed  those  for  the  vibrational  transitions  of  hydrogen  or  nitro¬ 
gen.  Hence  formula  (11-52)  gives  an  upper  limit  to  the  recombination  coef¬ 
ficient  provided  the  water  vapor  has  not  such  low  abundance  relative  to 
the  main  gas  that  it  does  not  control  the  conductances  on  the  outer  side 
of  the  minimum.  This  upper  limit  is  of  interest  in  connection  with  the 
problem  of  recombination  to  complex  ions  which  are  such  that  the  energy 
levels  of  X  lie  much  closer  together  than  the  hydrogenic  levels  assumed 
in  the  semi-quantal  calculations. 

An  experimental  study  of  the  decay  of  Pb+  ions  in  flames  at 
atmospheric  pressure  has  been  carried  out  by  Hayhurst  and  Sugden  (1965). 
Mass  spectroscopic  measurements  showed  that  the  ratio  [Pb+]/[Pb~ . HtO]  was 
about  300  which  enabled  Hayhurst  and  Sugden  to  dismiss  the  dissociative 
recombination  process 

Pb*.Ht0  +  e->Pb  +  H,0  (11-53) 

and  to  infer  that  the  decay  was  caused  by  recombination  following  a  col  - 
lision  between  an  electron  and  a  neutral  molecule,  probably  ^0  as  in 
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Pb+ + e  +  H  a0  ->■  Pb + HjO. 


(11-54) 


Their  results,  expressed  as  effective  binary  recombination  coefficients, 
are  given  in  Table  11-11  where  they  are  compared  with  the  corresponding 
result  obtained  by  the  semr-quantal  method  and  by  the  classical  diffusion 
method.  It  may  be  seen  that  the  measured  values  of  the  recombination 
coefficient  are  straddled  by  the  calculated  values  in  the  predicted  manner. 
In  view  of  the  ion  concerned  it  is  reasonable  that  they  should  lie  closer 
tc  the  upper  limit  than  to  the  lower  limit. 

“I-  4* 

Kelly  and  Padiey  (1969)  have  investigated  rhe  decay  of  Ga  ,  In 
and  Tl+  ions  in  flames  at  atmospheric  pressure.  In  accordance  with  the 
theory  the  recombination  coefficients  obtained  for  these  heavy  ions  did 
not  differ  significantly  from  each  other  or  from  Lhe  recombination  coef¬ 
ficients  for  Pb+  ions  (Table  IT- 1 1 ) - 

Unfortunately  direct  determinations  have  not  been  made  of  the 
recombination  coefficients  for  Na+  or  other  light  ions  the  parent  atoms 
of  which  are  nearly  hydrogenic.  Hu.iever,  there  have  been  studies  (Hay- 
hurst  and  Sugden  1965,  1967;  Kelly  and  Padiey  1969)  of  the  growth  of  the 
ionization  of  alkali  atoms  in  flames  at  atmospheric  pressure.  The  growth 
is  attributed  to  processes  like 


TABLE  11-11.  Recombination  coefficients  for  Pb4-  ions  in  damp 

NITROGEN-HYDROGEN  MIXTURES  AT  ATMOSPHERIC  PRESSURE 

recombination  coefficient  cm1  s  ' 1 


composition  of  mixture! 
volume  percentage 


calculated 


temperature, 

, - 

-  ' - 

, - 

~V - 

t?,K 

X. 

H. 

H.t! 

measured! 

lower  limit 

npj>f»r  limit 

1CK0 

67 

14 

Ml 

2.1  - 

9.2  * 

*2/2  * 

1785 

63 

iti 

2i 

1.2  ’ 

s.t>-> 

1.7*" 

1800 

63 

16 

21 

1.3  * 

7.9  * 

i.«;  * 

2005 

5V 

25 

s.l  » 

4.1  • 

1.1** 

2010 

50 

25 

6.2  * 

4.0-’ 

l.!“* 

2260 

47 

23 

Mt 

6.1  ’ 

3.1“* 

2270 

47 

23 

30 

6  1  * 

3.0  * 

i»  S“* 

+  The  temperature,  the  recombinat.on  coefficient  anil  the  percentage  of  water  vnjKiur  arc 
given  in  Hayhurst  &  Sugden  (1965);  other  information  on  the  composition  is  gnen  in  Hav hurst 
A  Sugden  (1966,  1967). 

The  indices  give  the  powers  of  It*  by  which  the  entries  should  i>o  multiplied. 


(11-55) 


Na + Ht0  ->  Na+ + e  +  HtO, 


and  from  the  measurements  carried  out  the  recombination  coefficient  for 
the  inverse  process 

Na+  +  e  +  HjO-*-Na+HjO  (11-56) 


may  be  deduced-  The  values  obtained  are  almost  the  same  as  those  for 
the  other  species  of  ion  mentioned  in  the  two  preceding  paragraphs.  In¬ 
deed  Kelly  and  Padley  (1969)  comment  that  there  is  no  detectable  dependence 
of  the  recombination  coefficient  on  any  physical  property  peculiar  to  each 
element  (such  as  the  ionization  potential  or  the  electronic  structure). 

The  comment  is  not  in  accord  with  the  expectation  that  the  recombination 
coefficient  should  lie  closer  to  the  lower  than  to  the  upper  limit  (see 
Table  II-H)  if  the  energy  levels  of  system  X  are  almost  hydrogenic.  This 
does  not  necessarily  mean  that  the  physical  picture  presented  here  is  in¬ 
correct;  it  may  merely  mean  that  the  semi-quantal  method  underestimates 
the  recombination  coefficient  for  a  hydrogenic  system  by  a  factor  of  about 
2  which  in  theory  would  not  be  surprising. 

Further  measurements  are  needed,  especially  on  recombination  in 
different  gas  mixtures. 
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III.  THE  SCATTERING  OF  THERMAL  ELECTRONS  BY  CARBON  MONOXIDE 
by  A.  Dalgarno 


ABSTRACT 


Close-coupling  calculations  are  carried  out  for  the  scattering  of 
electrons  with  energies  up  to  0.1  eV  by  a  rigid  rotator  with  dipole  and 
qoadrupole  moments  equal  to  those  of  carbon  monoxide.  By  appropriate 
choice  of  the  short  range  potential,  the  experimental  data  on  CO  momentum 
transfer  cross  sections  can  be  closely  reproduced.  For  the  scattering  of 
thermal  electrons,  CO  is  similar  to  Ng  with  the  addition  of  a  permanent 
dipole  moment.  The  observed  cross  section  minimum  occurs  through  the 
addition  of  the  decreasing  electron-dipole  contribution  and  an  increas¬ 
ing  spherically  symmetric  short  range  contribution. 

Detailed  angular  distributions  are  presented  for  elastic  and  rota¬ 
tional  excitation  collisions.  Transitions  in  which  the  rotational  quantum 
number  changes  by  unity  are  strongly  peaked  in  the  forward  direction  and 
transitions  in  which  the  rotational  quantum  number  changes  by  two  are 
nearly  isotropic.  The  elastic  scattering  is  peaked  in  the  backward  dir¬ 
ection. 
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A.  INTRODUCTION 

The  cross  section  for  inelastic  scattering  of  low  energy  electrons 
in  polar  gases  resulting  in  excitation  of  rotational  modes  of  the  polar 
gas  is  large,  for  polar  molecules  with  sufficient  static  dipole  moment. 
Knowledge  of  the  cross  section  is  required  in  order  to  study  the  energy 
loss  process  of  atmospheric  electrons,  and  in  order  to  determine  the 
rotational  temperature  of  the  atmospheric  polar  molecules.  Explicit 
calculations  have  been  carried  out  for  low  energy  electrons  in  CO.  Momen¬ 
tum  transfer  and  rotational  excitation  cross  sections  have  been  obtained, 
including  angular  differential  cross  sections. 

B.  BACKGROUND 

There  has  been  much  discussion  of  the  scattering  of  thermal  electrons 
by  polar  molecules  (see  Christophurou  and  Christodoulides  1969,  Christo- 
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phorou  and  Pittman  1970).  A  case  of  particular  interest  is  that  of  carbon 
monoxide  because  swarm  data  are  available  that  have  yielded  scattering 
cross  sections  over  a  considerable  range  of  electron  energies  (Hake  and 
Phelps  1967).  An  accurate  theoretical  description  is  afforded  by  the 
close-coupling  formalism  of  Arthurs  and  Palgarno  (I960),  previous  calcu¬ 
lations  of  the  momentum  transfer  and  rotational  excitation  cross  sections 
for  thermal  electrons  in  CO  having  been  restricted  to  the  first  Born 
approximation  (Singh  1970).  Close  coupling  calculations  for  simultaneous 
rotational  and  vibrational  excitation  at  energies  above  0.4  eV  have  been 
reported  (Itakawa  and  Takayanagi  1969) . 

C.  THEORY 

We  assume  that  for  the  scattering  of  thermal  electrons  the  molecule 

can  be  regarded  as  a  rigid  rotator  specified  by  rotational  and  projection 

quantum  numbers  j  and  nr .  In  the  close-coupling  formalism  the  angular 
momentum  j  of  the  rotator  is  compounded  with  the  orbital  angular  momen¬ 
tum  1  of  the  electron  to  form  a  total  angular  momentum  J.  The  interaction 
potential  V(r),  where  r  is  the  position  vector  of  the  electron  measured 

from  the  centre  of  mass  of  the  rotator  (see  Figure  III-l),  has  matrix  elements 

<jT; J\V\j’l’:Jy,  diagonal  in  J.  between  the  different  channel  angular 
momentum  eigenfunctions. 

The  scattering  matrix  SJ(jl',j'l')  is  defined  by  the  requirement  that 
the  solution  uf-/r(r)  of  the  set  of  coupled  equations 


*2  d3  /'(/'  + 1)  ,2]  Jli(  , 


2 fi  l  dr3 


4.  y  v  j  r,j\v\fr;J>u,-r{r)  =  o 


(in-i) 


where  kf-f  is  the  channel  wave  number,  behaves  asymptotically  as 


uJJ Ur) -8 exp{  -  i(knr  —  A/nj) 

“  (j^J  S'V1'  >T>  exPt  +  '(A,  /  - i M} 


(Arthurs  and  Dalgarno  1960) . 


Then  the  differential  cross  section  for  scattering  through  an  angle  0 
with  excitation  from  level  j  to  level  j7  is 


<-!)'-'  C- . „  ,  » 

— df — 


(III-2) 


where  px  are  Legendre  polynomials  and  Ax  are  expressible  in  terms  of 
the  scattering  T  matrices 

TJ(jhjl')  =  8„.8ir-S’(ji;j'n- 


Figure  J1I - 1 .  The  coordinate*  r  and  S'  are  shown.  M  i*  the  centre  of  ma**  of  CO. 


Introduce  the  algebraic  coefficient 


Z(abcd;  ef)  = 

/a  c  f\  (a  b  e  \ 
■  lW2r+  11  (2d+  Y\(2f  +  llll,a  I  J  \  ! 


{(2a+l)(2A+  l)(2r+ 1) (2d+  1)(2/+  l)}l,a  ‘  ^  *J 

,  (a  c  f\  la  b  e) 

where  0  0 j  and  (</  c  f  \  are  resPectively  3-j  and  6-j  coefficients. 


=  i  £  Y  Y  Y  Y 

j,-0  j,-o  I..I/.-/I  I.-I X-r.  ir-i/,-1-  V«l7,~n 

x  Z(l1J1l3J2;jX)Z(l1'J1la'Ja-  j'X) 

*T,'U'W\jhYTHrt3'-,jl2). 


The  total  cross  section  for  the  j-^j'  transition  is  given  by 


(III-3) 


42 


B  i 


R 


E 


•lR  V^O&fjyjjvg-^-,; 


-*•/> 


<2i+  DV 


«  /+*  /+r 


2  2  2  (2/+l)!r^V';;7)l2 


(2/+  l)^//a /-o  i-ii-ii  r-u-n 


and  the  momentum  transfer  cross  section  is  given  by 


•Jti>-  2°n(J  -+)')  =  J  ^^^(l-COS^df 


(2/+l)V7 


2^oO->/)-MiOW')}- 


D .  INTERACTION  iOTENTIAL 


To  construct  an  interaction  potential,  we  assumed  initially  that 

•18 


apart  from  its  permanent  electric  dipole  of  moment  0.112  x  10  esu  cm, 
carbon  monoxide  is  essentially  similar  to  the  isoelectronic  nitrogen 
molecule  (Hake  and  Phelps  1967).  The  interaction  potential  can  be  written 
as  a  sum  of  a  polarization  interaction  Fp(r),  a  dipole  interaction  Va(r) 
and  a  quadrupole  interaction  .  For  Vp(r)  we  adopted  the  form 


VJlr)  =  - 


oo*2 


r  ^  r. 


2(r*  +  rc*)* 

(r-r»)3 


2 (r*+r*)*  2(r»  +  r0a)26p2  +  (r-rp)3 


P2  (cos6‘) 


r  £  rp 


where 


*o  =  $(z,  +  2s:J 


«2  = 


a,  and  being  respectively  the  parallel  and  perpendicular  polariza¬ 
bilities  (see  Hirschfelder  et  al.  1954),  and  r0,  rp  and  bv  are  disposable 
parameters.  The  argument  r  is  the  position  vector  of  the  electron  rela¬ 
tive  to  the  center  of  mass  of  CO,  and  0 7  is  the  angle  between  r  and  the 
molecular  axis.  For  Va(r)  we  adopted  the  form 
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Fd(r)  =  0  r  <  rd 

Fd(f)  =  --  (r~rA — _  p  (cos#')  r  2  ra 
r*bS+(r-rJ* 

where  D  is  the  dipole  moment  of  carbon  monoxide,  and  for  Va(*)  the  form 


F„(r)  =  0  r  <  r, 
y  (r\  (r  ~  rq)2-—  P2  (co *B‘) 

V*T)  r3  b*  +  (r-r*Y 


where  Q  is  the  quadrupole  moment  of  carbon  monoxide  (the  <?«  of  Stogryn 
and  Stogryn  1966),  and  ra>  ba,  rQ  and  bq  are  disposale  parameters. 

The  potentials.  Fp,  Fd  and  F,  have  the  correct  asymptotic  behaviors 


^p(r)~-i^-i^Pa(cosO') 

VM-^VAcosO') 

r3 

F„(r)~  ^Pa(cosr). 

The  form  of  the  spherically  symmetric  part  of  Vp{r)  is  that  which  proved 
successful  in  the  interpretation  of  the  scattering  of  low  energy  electrons 
by  molecular  nitrogen  (Sampson  and  Mjolsness  1965)  and  its  short  range 
behavior  reflects  in  a  phenomenological  way  the  effect  of  exchange  forces. 
The  elastic  sc  -.ctering  cross  sections  depend  initially  on  the  short  range 
spherically  symrne  .ric  component  of  F(r)  and  we  chose  r0  so  that  our  calcu¬ 
lations  reproduced  the  momentum  transfer  cross  section  measured  at  the 
energy  of  0.03  eV. 

The  interaction  jf  an  electron  with  carbon  monoxide  cannot  be  properly 
described  by  a  local  velocity  independent  potential  for  values  of  r  smaller 
than  the  equilibrium  C-0  separation  for  r,  .  The  cross  sections  for  thermal 
electrons  are  however,  insensitive  to  the  detailed  form  of  the  orientation- 
dependent  terms  for  r^r,  and  we  imposed  arbitrarily  the  requirements 
that  Fd  and  F,  vanish  for  r  $  ir.and  r  <  \r,  respectively,  these  values  being 
our  estimates  of  where  the  corresponding  components  of  the  static  field  of 
CO  change  sign.  We  terminated  the  angular  term  in  Fp  also  at  r  =  |r,  . 
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The  values  adopted  for  the  potential  parameters  are  listed  in  Table  III-I. 
The  dipole  polarizabilities  are  taken  from  Bridge  and  Buckingham  (1966), 
the  dipole  moment  from  Barnes  (1958)  and  the  quadrupole  moment  from 
Stogryn  and  Stogryn  (1966). 

E.  CALCULATIONS 

Because  of  the  long  range  of  the  electron-dipole  interaction,  many 
partial  waves  contribute  to  the  scattering.  Even  close  to  threshold, 
over  100  values  of  J  must  be  retained.  The  computational  difficulty  can 
be  avoided  by  recognizing  that  as  J  increases,  the  exact  T  matrices 
rapidly  approach  the  first  Born  approximation  T  matrices.  Indeed  in  no 

TABLE  III-l .  value*  of  the  potential  parameters  in  V(r) 


Parameter 

Value 

*0 

1-977  A3 

*3 

0-355  A5 

0-693  A 

rP 

0-564  A 

b „ 

0-1  A 

D 

0-112x10  '*  esu 

T 4 

0-564  A 

ba 

0-1  A 

Q 

-2-5  x  10  ’*  esu  i 

N 

0-846  A 

b. 

0-1  A 

case  did  we  find  it  necessary  to  solve  the  close-coupled  equations  (III-l) 
for  J  >  7.  The  TJ  matrices  for  the  first  Born  approximation,  T  can 
be  easily  calculated.  With  the  long  range  electron-dipole  interaction, 
the  first  Born  approximation  allows  only  transitions  in  which  j  changes 
by  unity,  and  with  the  long  range  electron-quadrupole  interaction,  the 
first  Born  approximation  allows  only  transitions  in  which  j  changes  by 
two  (for  E  states  of  diatomic  molecules).  For  both  interactions,  the 
first  Born  dif fe.-.e-itiai  scattering  cross  sections  take  simple  forms  (see 
Crawford  et  al.  I  f.  ;) . 


The  values  of  Ax  for  large  values  of  X  are  needed  only  for  the  con¬ 
struction  of  the  differential  scattering  cross  section  at  small  angles. 
But  for  large  values  of  X,  Ax  and  Ax'  approach  the  Born  approximation 
.<4/(0)  .  Thus  if  /  and  F  are  sufficiently  large. 


do 

df 


O'-i'M 


dt  4(2/ +  l)fey>a  A-0 


where  doB/df  is  the  Born  approximation  to  the  differential  scattering 
cross  section. 


F .  RESULTS 

It  is  clear  from  the  comparison  of  the  first  Born  approximation  and 
the  fixed  rotator  approximation  of  Altshuler  (1957)  that,  except  very 
close  to  the  excitation  thresholds,  the  momentum  transfer  cross  section 
for  electrons  scattered  by  CO  is  nearly  independent  of  the  initial 
rotational  state  of  CO  (Crawford  et  al.  1967,  Bottcher  1970,  Crawford  1970) 
and  accordingly  of  the  gas  temperature.  Thus  the  close-coupling  results, 
obtained  from  the  numerical  integration  of  (III-l),  yielded  a  value  of 

7.41  X  ^  for  cr  at  0.1  eV  both  for  an  initial  j  of  0  and  of  10.  At 

3  ®  q2 

5  x  10  J  eV,  the  values  of  a  for  j  =  0  and  j  =  3  were  respectively  13. 5A 
2  m 

and  13.0  X  .  The  remainder  of  our  calculations  were  carried  out  with 
j  =  0. 

The  value  of  r0  that  reproduces  the  experimental  cross  section  at 
0.03  eV  is  0.693X.  Figure  III-2  compares  the  close-coupling  results  for  cr^ 
with  the  measured  values  at  other  energies  and  demonstrates  that  the 
harmony  between  theory  and  measurement  is  maintained.  Figure  III-2  also  in¬ 
cludes  the  Altshuler  approximation  (1957),  which  is  essentially  equivalent 
to  the  first  Born  approxkmation  corresponding  to  the  electron-dipole 
interaction  (Crawford  et  al.  1967).  At  low  energies  the  collision  is 
dominated  by  rotational  transitions  due  to  the  electron-dipole  interaction 
and  because  the  scattering  is  strongly  peaked  in  the  forward  direction  and 
the  dipole  moment  is  not  large,  the  scattering  is  accurately  described  by 
the  first  Born  approximation.  The  cross  section  due  to  the  electron¬ 
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Figure  III-2.  The  momentum  transfer  cross  section  for  electrons  in  CO.  Full 
curve  derived  from  swarm  data  by  Hake  and  Phelps  (1967);  broken  curve,  the 
Bom  approximation  of  Altshuler  (1957);  x,  the  results  of  our  close-coupling 

calculations. 


dipole  interaction  decreases  as  E  .  The  observed  minimum  in  cr  near 

m 

0.03  eV  occurs  because  at  higher  energies,  the  elastic  scattering  by  the 

spherical  component  of  V(r)  provides  the  dominant  contribution  to  a  . 

m 

The  contribution  from  quadrupole  excitation  to  j 1  -  2  levels  is  nearly 
independent  of  energy  and  it  is  not  negligible.  Thus  scattering  by  CO 
is  qualitatively  similar  to  scattering  by  with  the  addition  of  electron- 
dipole  scattering.  A  detailed  analysis  is  given  in  Table  III-2  which  lists 
the  individual  contributions  to  cr^  from  elastic  scattering  Aj  =  j'-j  =  0, 


<* 

TABLE  III-2  .  Momentum  transfer  cross  sections  in  A* 


Energy  (eV) 

0-005 

0-01 

003 

005 

0-1 

o.(0  -0) 

051 

0-97 

2-39 

3-51 

5-75 

o.(0  -*1) 

1 1  -83 

5 -Sc» 

1-90 

M2 

0-55 

o.(0  -*2) 

0-96 

103 

1-06 

1-07 

111) 

o.(0) 

13-30 

7-8f> 

5-35 

5-70 

7-40 

single  quantum  excitation  Aj  =  1  and  two  quantum  excitation  Aj  =  2.  There 
occur  some  contributions  to  Aj  =  0  and  Aj  =  2  transitions  from  the  virtual 
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dipole-dipole  sequences  but  for  CO  they  are  not  large  and  almost  all  of 
am<0  -*  0)  arises  from  l7P(f)  and  almost  all  of  cr^CO  -♦  2)  from  V^r)  . 

For  scattering  by  a  molecule  with  a  large  dipole  moment  such  as  CN  (Iti- 
kawa  and  Takayanagi  1969,  Crawford,  Allison  and  Dalgarno  1969),  the  situa¬ 
tion  is  reversed. 

The  elastic  and  excitation  cross  sections  cr(0~*j/)  are  listed  in  Table  III- 
3.  Hie  values  of  tr(O-l)  are  close  to  those  given  by  the  first  Born 
approximation. 


TABLE  III-3.  Total  elastic  and  excitation  I 


1  sections  in  A* 


Energy  (eV' 

0-005 

0-01 

0-03 

0-05 

0-1 

<7(0-0) 

0-40 

0-74 

1-78 

2-58 

4-15 

o(0-l) 

45-06 

26-87 

11-20 

7-34 

4-09 

o(0-2) 

0-96 

1-03 

1-07 

1-08 

1-11 

o(0) 

46-42 

28-64 

14-05 

11-00 

9-35 

The  differences  between  cr(j-*j')  and  ffjj-j')  are  a  consequence  of  the 
angular  distributions  do(j—j')ldr.  The  calculated  distributions  arc  illustrated 
m  Figures  III-3,  III-4  and  III-5.  Those  for  do(0  -* l)/df  are  given  in  the  form 
of  the  ratio  to  the  distribution, doB/df,  obtained  from  the  first  Born 
approximation  for  the  electron-dipole  interaction.  The  angular  distribu¬ 
tions  are  required  in  the  theory  of  recombination  in  molecular  gases 
developed  by  Bates  et  al.  (1970). 


0 

1  04; 
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Figure  III-3.  The  angular  distribution  of  electrons  of  various  energies  scattered  by 
CO  in  an  elastic  0—0  rotational  transition. 


100, 


S 

Figure  III-4.  Thc  angular  distribution  of  electrons  of  various  energies  scattered  by 
CO  in  an  inelastic  0-»  1  rotational  transition.  The  distribution  is  given  as  the  ratio 
to  the  angular  distribution  calculated  from  the  Bom  approximation  assuming  an 
ele'.tron-dipole  interaction. 


The  differential  cross  section  for  Aj  =  1  scattering  is  charply  peaked 
at  3mall  angles  behaving  there  almost  as  (l-kir  cosBIkj,)'1  ,  as  expected  for 
an  electron-dipole  interaction  (see  Takayacagi  1966,  Crawford  et  al.  1967), 
and  the  differential  cross  section  for  Aj  =  2  scattering  is  nearly  iso¬ 
tropic  as  expected  for  an  electron-quadrupole  interaction  (Gerjuoy  and 
Stein  1955).  The  elastic  scattering  cross  section  depends  upon  the  detailed 
form  of  V(r)  and  with  our  choice  it  peaks  in  the  backward  direction.  The 
(l-cos9)  factor  in  c suppresses  the  strong  forward  peak  associated  with 

Aj  =  1  scattering  and  cr  (0-1)  is  much  smaller  than  ct(0  -*  1).  The  cross 

m 

sections  are  comparable  both  for  Aj  =  0  and  Aj  =  2  transitions. 
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IV.  ASSOCIATIVE  DETACHMENT  OF  0  AND  0 
by  A.  Dalgarno  and  J.  C.  Browne 

Theoretical  studies  of  the  rate  coefficient  for  the  associative  de¬ 
tachment  reaction 

X  +  Y“  -  XY  +  e 

requires  accurate  knowledge  of  the  potential  energy  curves  for  both  the 
neutral  and  negative  ion  molecule.  Also  required  are  calculations  of 
the  half-width  for  the  radiationless  transition.  The  calculation  of  the 
half-width  is  very  difficult  for  complex  systems,  so  that  entirely  pre¬ 
dictive  calculations  have  been  carried  out  only  for  the  simple  reac¬ 
tion  (Ref.  1-3) 


H  +  H 


H2  +  e 


A  critical  parameter  in  the  calculation  of  the  rate  coefficient  for 
associative  detachment  of  0  and  0  : 

0  +  O’  -  02  +  e 

are  the  interaction  energies  at  the  intersections  of  the  C>2"  and  02  poten¬ 
tial  energy  curves.  If  any  one  of  the  crossing  points  lies  energetically 
below  the  asymptotic  limit  at  infinite  separation  ,  the  process  can  be 
expected  to  be  rapid  at  thermal  velocities. 


Accordingly  we  have  carried  out  elaborate  configuration  interaction 

calculations  for  the  ^1  and 

8 

previous  work. 


-u  status  of  0^  as  an  extension  of  our 


The  orbital  exponents  for  the  ?p  orbitals  were  chosen  so  that  the 
resulting  curve  had  an  electron  affinity  of  0.43  eV  with  respect  to  the 
Schaefesr-H arris  calculations  for  02  (J.  Chem.  Phys.  48,  4946,  1968),  when 
it  is  adjusted  so  that  the  difference  of  the  separated  atom  limit  and 
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the  Schaefer-Harris  07  curve  has  the  experimental  value.  The  0?  (  n  ) 

3  .  ^  ^  8 

crosses  the  0_(X  £  )  curve  at  a  sepavation  of  about  2.15  a0  at  an  energy 

“  g 

below  the  asymptotic  limit.  The  precise  value  of  the  :rossing  is  uncertain 

because  we  have  discovered  shift  errors  in  the  published  Shaefer-Harris 

0,  curves  that  will  have  to  be  corrected.  However,  there  seems  to  be 
1  2 

no  douot  that  approach  along  the  II  curve  leads  to  associative  detach- 

g 

ment  at  thermal  velocities. 
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V.  PSEUDO-POTENTIAL  CALCULATION  OF  ATOMIC  INTERACTIONS 
by  A.  Dalgarno  and  G.  A.  Victor 

ABSTRACT 

The  use  of  atomic  psuedo-potentials  in  calculations  of  atomic  inter¬ 
actions  is  briefly  discussed.  It  is  shown  that  an  appropriate  choice 
of  pseudo-potential  leads  to  cne  appearance  of  the  van  der  Waals  inter¬ 
actions  in  first  order.  Care  is  needed  in  the  choice  of  the  core-core 
interaction  potential  lest  spurious  long-range  forces  appear.  Explicit 
calculations  are  reported  for  the  aad  ground  states  of  L^+ 

and  comparison  made  with  the  results  of  conventional  calculations.  High 
accuracy  is  obtained  over  a  wide  range  of  internuclear  separations. 

A.  INTRODUCTION 
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Rydberg-Kel in-Rees  (RKR)  methods  yield  accurate  potential  energy 
curves  for  diatomic  molecules  for  spectroscopically  measured  attractive 
potential  curves  in  the  neighborhood  of  the  equilibrium  internuclear 
separation.  Similarly,  at  very  large  internuclear  separations,  ab-initio 
or  semi-empirical  methods  lead  to  detailed  knowledge  of  the  van  der 
Waals  terms  for  various  molecules.  For  intermediate  internuclear  sep¬ 
arations,  much  less  data  are  available.  Semi-empirical,  effective  poten¬ 
tial  methods  can  provide  accurate  data  for  selected  systems  at  intermed¬ 
iate  internuclear  separations,  where  other  methods  are  inappropriate. 

For  many  processes  of  atmospheric  interest,  such  data  are  a  prerequisite 
for  quantitative  predictions.  Preliminary  studies  of  the  molecular  poten¬ 
tial  curves  for  Li2+  have  been  carried  out  using  accurate  semi-empirical 
effective  potential  techniques.  The  results  are  very  encouraging,  as  the 
current  results  compare  very  favorably  with  very  elaborate  ab-initio 
calculations.  Such  procedures  can  be,  in  principle,  extended  to  more 
complex  systems  as  long  as  the  number  of  essential  valance  electrons  is 
not  large.  A  review  of  these  studies  is  given  in  the  following  discus¬ 
sions.  Further  preliminary  calculations,  for  alkali  systems  on  inert 
gases  have  been  carried  out  which  demonstrate  the  utility  of  the  methods 
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for  scattering  problems.  Using  these  methods,  limited  to  one  valance 
electron,  studies  can  be  made  of  charge  transfer  processes  in  alkali- 
alkali  ion  systems  as  well  as  quenching  processes  in  alkali- inert  gas 
systems.  Treatment  of  up  to  two  valance  electrons  is  very  possible, 
leading  to  the  possibility  of  studying  many  diverse  atomic  scattering 
problems;  however,  it  seems  that  relevant  data  for  few  processes  of 
interest  for  normal  and  disturbed  atmospheres  can  be  obtained  without 
encountering  severe  difficulties. 

B.  DETAILED  DISCUSSION 

Psuedo-potential  or  model  potential  methods  have  been  successfully 
applied  to  several  problems  of  quantum  chemistry  (Ref.  V-l)  and  they  can 
also  be  of  value  in  the  determination  of  interatomic  forces  at  intermediate 
and  large  nn-l.ear  separations.  Consider  a  Rydberg  electron  with  quantum 
numbers  n£  and  position  vector  ra  moving  in  the  field  of  a  spherical  core 
with  excess  charge  za  located  at  the  origin  A  of  coordinates  and  construct 
an  effective  single  channel  potential  Va(ra)  with  eigenvalues  equal  in  mag¬ 
nitude  to  the  observed  ionization  potentials B^{nl) of  the  Rydberg  states. 
Suppose  that  this  atomic  system  is  interacting  with  another  spherical 
core  B  at  a  distance  R  from  A  and  that  the  binding  energies  of  the  val¬ 
ance  states  associated  with  the  second  core  can  be  reproduced  by  an 
effective  potential  V'b^b)  •  The  problem  of  calculating  the  adiabatic 
interatomic  potential  t(R)  reduces  to  finding  the  appropriate  solution  of 
the  one-electron  equation  (written  in  atomic  units) 


{-  Va(ra)  +  V’b(rb)+  ra,  rb)-<(R)l  *(ra)  =  0  . 


(V-l) 


where  Vc  represents  the  effects  of  the  core-core  interaction.  At  large 


distances  R  ,  the  system  dissociates  to  a  Rydberg  state  of  either  A  or 
B. 


Suppose  we  are  interested  in  the  state  that  dissociates  into  the 
n2  state  of  A  with  wavefunction  «a(n/  ra).  such  that 
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{  -  ;.  v2  +  Va{ra)  -  £a(nO}  ^(k/)  =  0  . 


(V-2) 


then  at  large  distances  R  ,  to  first  order. 


e(R)  =  (<PA{nl  ira)iVb(rb)+Vc(^,  ra,  rb)|  <1>a(nl\ra.):  . 


(V-3) 


If  we  include  a  polarization  term  in  the  core  potential,  then  asymptotically 


Z  b  «b 

Vb<rb>  v 

*T>  2rg 


(V-4) 


where  is  the  dipole  polarizability  of  the  core.  Hence  for  spherically 


symmetric  systems. 


<(R)  =•  \  is  V  (R,  r  r  )  ns  .  £.  ^  -  j?  (ns  r2  ns  . 

C  a  b  R  2r4  «6  a 


(V-5) 


It  is  cle»  ‘"’at  the  core-core  potential  Vc(#.  ra.rb)  must  contain  asymp¬ 
totically  ’  only  the  nuclear  term  ZaZb.  R  but  also  the  polarization 
term  (--•  »'  <yb  R*  .  The  omission  of  the  polarization  term  would  be 

particularly  serious  for  a  neutral  core  with  Z,  =  0.  A  olausible  choice 

b 

for  Vq{R,  ra.  r  b)  is 


(Z'a - 1>  a^b *  (2a -  n2  Vb(R>-  Vb(fii 


(V-6) 


if  the  molecular  system  idssociates  to  a  Rydberg  state  of  A. 

The  R  ^  term  in  e(R)  is  the  van  der  Waals  interaction.  It  here 
appears  in  first  order,  rather  than  in  second  order  as  in  conventional 
approaches,  and  it  is  easily  reproduced  by  variational  solutions  of  the 
one-electron  eigenvalue  equation  for  e(R).  The  van  der  Waals  interaction 
can  be  written  in  the  firm 


(V-7) 


2R^ 


ffnft 


where  the  r s  are  oscillator  strengths  and  the  e's  are  atomic  eigenvalues. 

If  we  assume  that  the  transition  <-.iergies  (e^  -  «  )  of  the  Rydb  rg  system 

A  are  much  smaller  than  those  (e^  “  e^)  of  the  core  B,  it  follows  that 

u  n 

(Ref.  V-2) . 


2  6 
C  ~  -  u-ins  r  ns  R 
D  a 


(V-8) 


Our  choice  of  pseudo-potential  leads  to  the  same  asymptotic  interaction. 

The  numerical  value  given  by  (V-8)  is  -3.39,  whereas  the  correct  value  is 
-3.32  (Ref.  V-3) . 

Conventional  calculations  of  the  interaction  energies  of  have 

been  reported  by  Bardsley  (Ref.V-4)  and  they  provide  a  valuable  quantitative 
test  of  our  model.  We  adopt  the  representation  in  atomic  units 


0.09615  . 


0.0073 


VL,-(r)  =  VHF(r)-  ^4  U  -  expf  -  (r  0.47)*)}  -  •  -  (l  -  exp[  -  (r/0.44)8]}  (V-9) 


where  vHF(r)  is  the  Hartree-Fock  core  potential.  This  potential  (V-9)  has  a 
2s  eigenvalue  of  -0.19824  au  compared  to  the  experimental  value  -0.19814  au 
and  a  2p  eigenvalue  of  -0.13039  aj  compared  to  the  experimental  value 
-0.13024  au.  An  accurate  2s  -  ap  energy  difference  is  important  in 
ensuring  that  the  interactions  nave  the  correct  (second-order)  asymptotic 
form  -1/2  a/R  where  a  is  the  dipole  polarizability  of  neutral  lithium. 

The  pseudo-potentiai  gives  a  2s  -  2p  energy  difference  of  0.06785  au 
compared  to  the  experimental  difference  0.06790  au. 

We  have  solved  the  one-electron  eigenvalue  equations  for  e(R)  for 
the  a  (2s)  and  >s^(2s)  orbitals  by  standard  variational  procedures  except 
that  the  molecular  integrals  involving  the  pseudo-potentials  were  eval¬ 
uated  by  numerical  quadrature.  The  results  are  compared  with  those  of 
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Bardsley  in  Figure  V-l.  The  agreement  is  close  over  a  wide  range  of  sep¬ 
arations.  For  the  ground  state  of  we  obtain  a  dissociation  energy 

of  1.30  eV  at  our  equilibrium  separation  of  6.0  oq  whereas  Bardsley  (1970) 
obtains  1.22  eV  at  5.8  ao-  Asymptotically  we  find 


«(ft)  »- 


80.7  2926 
ft4  '  ft6 


(V-10) 


The  choice  of  l'c(R-  ra.  rb)  is  not  unique.  Indeed  of  particular 
importance  for  1*0  ,  a  more  accurate  representation,  correct  at  large 
R.  is  provided  by 


>  = 


R 


-A  V  (R  >  -  Z.V  ( R 1- 
a  b  b  a 


^b’va  ^  *  f a 


a >  (r  )- 
a  a 


•a,yb 


R  r 


»K3 

R  rb 


b  ,  .  ^a^b  2b°a 

'*’b  ^  2R4  2  ft4 


(V-ll) 


t 
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1 


1 

§ 

2 


-J- 

where  a»a(r).  Jb(r)  are  suitable  cut-off  functions.  In  the  case  of  Li  we 
took 


w(r)  i  1  .  exp  (  -  (r  0.47)6|  . 


(V-12) 


as  suggested  by  Equation  (V-9).  Asymptotically,  c (R)  behaves  correctly 
according  to 


t{R)>  {nl  r2,l->  P0[cos6  ))‘nl 
no  a  *  t  a 


(V-13) 


(see  Baylis,  Ref.V-5).  For  o  the  P„  term  vanishes.  We  have  employed 

1  2  + 

the  more  complicated  form  of  V  in  a  calculation  of  ( (R)  for  the  £ 

C  g 


state  at  R  =  6.0  <j0 
i .  28  eV . 


The  binding  energy  is  changed  from  1.30  eV  to 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

Detailed  studies  of  atmospheric  behavior  following  nuclear  or 
natural  (aurora  or  PCE  events)  disturbances  requires  accurate  data  on 
cross  sections  and  rate  coefficients  for  a  highly  diverse  ensemble  of 
atomic,  molecular,  chemical,  and  transport  processes.  In  recent  years, 
experimental  studies  and  in  a  more  modest  way,  theoretical  investiga¬ 
tions  supported  by  DNA  have  provided  essential  input  data  for  subsequent 
computer  modeling  of  atmospheric  response.  Owing  to  the  relatively  wide 
range  of  possible  initial  conditions,  the  complexity  of  the  overall 
problem,  and  the  inevitable  delay  between  identification  of  important 
processes  and  the  subsequent  generation  of  accurate  input  data,  under¬ 
standable  time  delays  occur  between  input  data, generation  data  and 
application  to  computer  coding.  As  a  consequence,  data  requirements  and 
thus  emphasis  on  specific  processes  change  periodically.  Depending  on 
the  accuracy  requirements  imposed  on  the  input  data,  decisions  on 
whether  to  employ  experimental,  theoretical,  or  a  combination  of  approaches 
become  important  from  both  time  and  economy  criteria.  Such  decisions  are 
complicated  by  the  fact  that  equilibrium  thermodynamic  conditions  are 
seldom  encountered,  so  that  detailed  cross  section  data  are  usually  re¬ 
quired. 

GCA  Technology  Division  has  for  several  years  been  performing 
theoretical  studies  directed  towards  the  satisfaction  of  these  require¬ 
ments.  Emphasis  has  been  directed  towards  the  study  of  those  processes 
which  control  the  degree  of  atmospheric  ionization,  particularly  for 
those  processes  which  are  important  for  recombination.  Basic  cross  sec¬ 
tions  have  been  investigated  rather  than  reaction  rates,  in  order  to  cope 
with  the  non-thermal  equilibrium  requirements  of  many  applications. 

A  major  portion  of  our  investigations  has  been  concerned  with  the 
process  of  three-body  ion-electron  recombination  in  dense  gases  at  low 
temperatures  (Section  II).  In  such  complex  systems,  one  must  resort  to 
semi-quar.tal ,  statistical  theories,  which  should  lead  to  results  of 
acceptable  accuracy  for  the  necessary  applications.  The  studies  per- 
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f-nned  to  date  yield  sufficient  data  for  useful  extensions,  without 
extensive  calculation,  for  predictions  for  various  ions,  temperatures, 
third  body  concentrations,  and  third  body  systems  expected  to  be  of 
importance.  Further  significant  progress  leading  to  increased  accuracy 
can  only  be  obtained  by  a  major  reformulation  of  the  theory  requiring 
significant  theoretical  progress  in  complex  many-body  problems.  Such 
progress  is  unlikely  in  the  near  future  so  that  additional  extensive 
theoretical  study  of  the  process  is  not  indicated  at  this  time.  The 
body  of  data  generated  to  date,  serves  as  a  useful  basis  for  interpo¬ 
lation  and  extrapolation  for  systems  not  yet  studied. 

The  results  presently  available  for  study  of  the  process  of  assoc¬ 
iative  detachment  are  preliminary  in  nature.  The  calculations  for  the 
potential  energy  curves  for  0+0  should  be  continued  in  order  to  pro¬ 
vide  mote  accurate  data,  especially  in  the  region  of  the  important  curve 
crossings.  Similar  studies  should  be  performed  for  the  N  +  0  system. 
Only  an  upper  limit  for  the  cross  section  for  associative  detachment  can 
be  obtained  from  the  potential  energy  curve  data.  Accurate  quantitative 
predictions  require  the  calculation  of  certain  coupling  matrix  elements. 
Two  methods  can  be  suggested  for  the  evaluation  of  these  necessary  matrix 
elements.  One  involves  the  evaluation  of  the  electronic  molecular  Hamil¬ 
tonian  for  complex  energies  whereas  the  other  involves  the  direct  eval¬ 
uation  of  the  coupling  matrix  elements  using  single  center  molecular 
eigenfunctions.  Both  methods  are  within  the  state-of-the-art  for  large 
scale  molecular  calculations,  and  should  be  investigated  in  greater  de¬ 
tail  for  possible  application  to  the  problem  of  associative  detachment. 

The  studies  of  the  use  of  effective  potential  methods  for  atomic  and 
molecular  processes  has  provided  highly  accurate  results.  However,  in 
the  present  form,  they  are  of  limited  use  for  important  atmospheric 
species  and  progresses,  so  that  more  extensive  studies  are  not  indicated 
for  future  progress. 

Two  areas  of  study  should  be  introduced  into  future  programs.  The 
first  involves  the  use  of  classical  and  semi-quantal  methods  ter  the 


theoretical  investigation  of  several  classes  of  ion-molecule  reactions. 
This  study  is  indicated  at  the  present  time  because  the  classical  and 
semi-quantal  methods  appear  to  yield  accurate  results  for  many  processes, 
and  ion-molecule  reactions  form  an  important  class  of  reactions  in  the 
overall  understanding  of  disturbed  atmospheres. 

Higher  energy  processes  are  becoming  of  increased  importance  in  dis¬ 
turbed  atmosphere  studies.  Because  of  recent  theoretical  progress, 
quantitative  studies  of  the  process  of  diaelectronic  recombination  are 
becoming  possible.  As  a  consequence,  a  renewed  study  of  diaelectroni.: 
recombination  should  be  nitiated. 
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cule,  very  large  recombination  rates  will  result.  The  calculations  show 
that  the  rate  for  three-body  ion-electron  recombination  is  very  sensi¬ 
tive  to  the  concentration  of  water  vapor. 

The  process  of  associative  detachment 

X~  +  Y  -  XY  +  e 

is  an  important  mechanism  for  the  removal  of  negative  atomic  and  molec¬ 
ular  ions.v  Only  sparse  accurate  experimental  data  exist,  and  quantita¬ 
tive  theoretical  data  exist  only  for  H  +  H  collisions.  Quantitative 
prediction  of  the  rate  coefficient  for  associative  detachment  requires 
accurate  potential  energy  curves  for  the  negative  molecular  ion  XY  and 
the  neutral  molecule  XY,  especially  in  any  regions  of  curve  crossings 
involving  the  lower  levels  of  the  molecules.  Also  required  are  certain 
non-radiative  coupling  matrix  elements  between  the  electronic  states  of 
the  molecules.  The  theoretical  calculation  of  the  necessary  accurate 
potential  energy  curves  and  matrix  elements  represents  a  relatively 
difficult  theoretical  problem  for  species  of  interest  in  the  atmosphere. 
From  knowledge  of  the  potential  energy  curves,  useful  upper  limits  for 
the  rate  of  associative  detachment  can  be  obtained.  Improved  calcu¬ 
lations  of  potential  energy  curves  0^,  relative  to  .theoretical  Oj  curves, 
have  been  performed, for  the  important  atmospheric  process 


0  +  0  —  O2  +  e 


the  cross  section  could  be  as  large  as  4.6  ra 

o 

The  cross  section  for  rotational  excitation  in  electron  polar  mol¬ 
ecule  collisions,  * 

e  +  M(J)  e  +  M  (J  +  -AJ)  , 

where  M(J)  is  a  polar  molecule  in  rotational  state  J,  is  large.  Rota¬ 
tional  excitation  processes  represent  an  important  energy  loss  mechanism 
for  the  thermalization  of  supra-thermal  electrons  produced  by  an  atmo- 


spheric  disturbance.  Since  most  recombination  processes  proceed  much 
more  rapidly  for  low  energy  electrons  than  for  supra-theraal  electrons, 
the  rotational  excitation  process  plays  an  important  role  in  the  rate 
of  reduction  of  ionization  in  a  disturbed  atmosphere.  The  cross  section 
for  this  process  constitutes  an  essential  input  data  point  for  the  studies 
of  three-body  ion-electron  recombination  in  polar  gases.  Detailed  cal¬ 
culations  have  been  performed  for  the  cross  section  for  rotational  exci¬ 
tation  of  carbon  monoxide  by  electron  impact  using  close-coupling  methods. 
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